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ABSTRACT 



The Global Positioning System (GPS) is revolutionizing the science of navigation. Never before has 
there been a system that could provide real time, world wide, continuous coverage with the such precision. 
Yet, the accuracy achievable with GPS alone is not sufficient to achieve autonomous flight of an Unmanned 
Aerial Vehicle. However, when integrated with an Inertial Navigation System and other non-inertial sensors 
using a Kalman Filter, GPS supplies the critical positioning information to permit such an achievement. 

This thesis presents the selection and integration of a GPS receiver using Differential GPS (DGPS) 
in support of a UAV autonomous flight project. Contemporary electronic navigation systems are surveyed, 
GPS operation is reviewed, and a Motorola PVT-6 GPS receiver selected. Using the Motorola Proprietary 
Binary Format protocol, several software drivers were written in C to interface the information to an Intel 
80486DX CPU using the RS-232 serial communication standard. Finally, an examination is made to 
determine the maximum reacquisition time, the DGPS accuracies achievable and the effects of pseudorange 
correction latency on DGPS accuracy. 
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INTRODUCTION 



The last of a 24 satellite constellation was launched into 
orbit in June 1993, and with it, "the science of navigation is 
poised to take a great leap forward" (Carey, 1993, p. Bl) . The 
Global Positioning System (GPS) revolutionized electronic 
navigation by providing worldwide, real-time, continuous 
coverage with three dimensional accuracies at least an order 
of magnitude better than any other system in existence. Even 
with these capabilities, it is not the panacea some make it 
out to be. Absolute accuracies (i.e., accuracies using a 
single receiver) better than ten meters within two standard 
deviations cannot be achieved. 

Differential GPS (DGPS) shows great promise by further 
improving accuracies to within one to three meters. DGPS uses 
two GPS receivers where one is placed in a known location and 
the other is allowed to move. By computing the difference 
between the stationary receiver's actual and GPS determined 
locations, and then transmitting these differences to the 
other GPS receiver, the other GPS receiver can correct its 
output to obtain the improved position information. This 
improvement would be sufficient for an aircraft to achieve an 
autonomous landing with GPS alone were it not for the 
unreliability of the GPS signals in high speed, critical 
maneuvering evolutions. It is not uncommon for a GPS receiver 
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to lose track of one or more of its satellites for a several 



seconds, and thereby, be unable to provide positioning 
information. While not a catastrophe for a slow moving surface 
vehicle, an outage of even a few seconds could prove fatal for 
an aircraft. 

In an attempt to achieve the accuracy and reliability 
required for autonomous flight in a navigation package light 
enough for a small unmanned aerial vehicle (UAV) , this 
research effort is part of a project where the GPS information 
is combined with the outputs of an Inertial Measurement Unit 
(IMU) and with data from other non-inert ial navigation sensors 
such as airspeed and altitude indicators in a Kalman Filter. 
A contribution is made by selecting the best GPS receiver, 
defining its method of employment and writing the necessary 
software drivers to smoothly interface them with a Central 
Processing Unit (CPU) . 

In the following chapter, contemporary electronic 
navigation systems are surveyed, discussing their advantages 
and disadvantages. Realizing the obvious benefits of GPS, 
Chapter III examines the User Segment of GPS emphasizing those 
details necessary to intelligently select and use the GPS 
receiver in a manner that best supports the various project 
requirements. Chapter IV brings it all together by detailing 
the receiver selection and capabilities as well as the 
criteria for its use. It also describes the hardware and 
software necessary to interface the receivers with the CPU. 
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Three tests were conducted to confirm expected operational 
capabilities. They are described and the results from them are 



presented in Chapter V. The appendices 
specifications, the ephemeris algorithm 
drivers written in C code as well as other 
the user and the programmer. 



include receiver 
and the software 
information to aid 
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II. SURVEY OF ELECTRONIC NAVIGATION METHODS 



The invention of wireless communications in the early 
twentieth century ushered in a new era in navigation allowing 
seafarers to safely pilot their ships in the absence of sight 
and sound for the first time in the history of mankind. The 
first use of radio waves required a vessel to transmit a 
continuous wave at some frequency. On shore, receivers with a 
rudimentary direction finding capability triangulated the 
ship's position and transmitted it back to the ship. The 
position was not very accurate by today's standards but served 
its purpose well in regions of the world where the sky was 
commonly obscured by clouds or fog when far from land. 

The next significant advance came with the invention of 
highly accurate chronometers. First, the quartz clock, and 
then, the atomic time standards with mind boggling stabilities 
providing accuracies better than one second in 30,000 years or 
one part in a trillion (10 12 ) (Logsdon, 1992, p.156). Such 
precision is required to allow the precise measurement of 
frequency. In addition, Transit (SatNav) and Global 
Positioning System (GPS) partially determine position by 
measuring the amount of time for the signal to travel from the 
satellite to the receiver. Radio waves travelling at the speed 
of light will advance one foot in one billionth of a second 
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(Kayton, 1990, p.336). It is not difficult to see the 

necessity for such accuracy in time keeping. 

In order to appreciate the quantum improvement in the 
science of navigation that GPS has achieved, it is worthwhile 
to review the major electronic navigation methods in use 
today. A brief review of some principles and characteristics 
of electromagnetic waves is given, followed by an individual 
description of each of the major navigation systems. While an 
effort has been made to succinctly describe the theory upon 
which the system operates, the important point is to notice 
the range, coverage and accuracy for each method. 

A. PRINCIPLES AND CHARACTERISTICS OF ELECTROMAGNETIC WAVES 

Radio waves in a vacuum theoretically travel at the speed 
of light. Wavelength is related to frequency by Equation (1): 

^ _ 2.997925 x 10 9 m/sec (1) 

where X is the wavelength in meters and / is the frequency in 
Hertz. Figure 1 demonstrates this relation for the 
electromagnetic spectrum. Information may be carried on a 
radio wave through one or more modifications to its shape. 
These modifications include changing the amplitude, frequency 
or phase of the radio wave and are appropriately called 
amplitude, frequency and phase modulation, respectively. In 
addition, a wave may be intermittently transmitted or pulsed. 
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This is called pulsed modulation. Figure 2 shows three 
representative examples of modulation techniques. 
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Figure 1. Wavelength Related To Frequency In The 
Electromagnetic Spectrum. (Hobbs, 1981, p. 217) 



Radio waves, like light, experience such phenomena as 
reflection, absorption, refraction, diffraction and 
interference. With respect to the earth's environment, the two 
most influential features on the distance and path radio waves 
travel are the earth's surface and the ionosphere. The 
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FREQUENCY MODULATED WaVE 




PULSE MODULATED WAVE 

Figure 2. Modulation Techniques Commonly Used. 

(Hobbs, 1981, p . 2 19 ) 

ionosphere has four layers and exists at a height of 
approximately 60 to 400 Km above the earth's surface. This 
highly charged and fluctuating segment of the atmosphere not 
only absorbs radio wave energy but also reflects it (Beck, 
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1971, p.95). Waves that are reflected by the ionosphere are 
called "sky waves". The ionosphere is affected by diurnal and 
seasonal fluctuations and by the level of sunspot activity 
which in turn affects the performance characteristics of sky 
waves. Sky wave phenomena is also greatly influenced by the 
frequency, angle of incidence, and the height and density of 
the various layers of the ionosphere. The sky wave frequencies 
most commonly observed are in the 300 KHz to 30 MHz band 
(Medium and High Frequency bands) . Because of the higher 
resistance of the earth's crust compared to the atmosphere, 
radio waves transmitted parallel to the earth's surface are 
slowed causing the wave to bend. The degree of bending is 
inversely proportional to the frequency. Thus lower 
frequencies will travel over the horizon in a direct path. 
These radio waves are called ground waves. As it might be 
expected, a receiver might travel along the surface of the 
earth from the transmitter origin and receive the ground wave, 
enter into an area shadowed by the earth's curvature with no 
signal, and then pick up the reflected "one-hop" sky wave even 
further out. Figure 3 illustrates this. Furthermore, when 
conditions are right, one-hop sky waves may experience 
reflection by the earth's surface and subsequent ionospheric 
reflection for a "two-hop" sky wave. Navigation systems such 
as Loran-C and Omega use reflected sky waves (when they exist) 
to extend their effective range. (Logsdon, 1992, pp. 95-98) 
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Figure 3. Sky Wave And Ground Wave Patterns. Ray 1 Is 
Refracted Through The Ionosphere, Ray 2 Is A One-Hop 
Skywave, Ray 3 Is A Two-Hop Skywave And Ray 4 Is A Ground 
Wave. (Hobbs, 1981, p. 225) 

Finally, several major navigation systems exploit radio 
wave phase measurement in order to determine two dimensional 
position. The principle of operation uses two synchronized 
transmitters located a fixed distance apart. The line between 
them is called the baseline. As a result, at half-wavelength 
intervals the phase difference between the two signals is 
zero. In this type of system, the half-wavelength segments are 
normally called lanes. As the distance from the baseline 
increases the lines of equal phase or "isophase lines" radiate 
outward from the foci of the two transmitters in hyperbolic 
fashion as shown in Figure 4. Thus, methods employing this 
technique are called hyperbolic radio navigation systems. If 
the two transmitters were to transmit a pulse instead of a 
continuous wave, then by measuring the time of arrival of each 
synchronized pulse, one's position might also be determined. 
The points where the times of arrival of transmitted pulses 
are the same describe a hyperbolic locus. This type of method 
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is also considered a hyperbolic radio navigation system. With 
the exception of radio beacons, all contemporary terrestrial 
radio navigation systems use some sort of hyperbolic 
navigation method. (Bowditch, 1984, pp. 1015-1016) 




Navigation Pattern. (Hobbs, 1981, p. 227) 

B. RADIO BEACONS 

Radio beacons are continuous wave transmitters modulated 
with a unique Morse code identification. They are usually low 
power transmissions in the low and mid frequency bands with a 
range of less than 200 miles. Most cannot be received beyond 
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20 miles (Hobbs, 1981, p. 231) . They operate using shipboard 
radio direction finding (RDF) equipment where two or more 
beacons are used to obtain a fix. The accuracy is highly 
dependent on the sensitivity of the RDF equipment as well as 
the range from the beacons. As Figure 5 illustrates, the 
geometry of the radio beacons relative to the ship may also 
play a significant role in the accuracy. Finally, the precise 
positions of the beacon transmitters must be known or the line 
of bearing is useless. In other words, the position 
information cannot be determined in terms of latitude or 
longitude until after it is plotted. 

C. LORAN-C 

Loran-C is the descendent of the first attempt at a 
hyperbolic navigation system developed during World War II for 
both ships and aircraft (Wilkes, 1987, pp. 31-34). From Figure 
6, it can be seen that it covers almost the entire northern 
hemisphere with only a very small area covered below the 
equator and then only by sky waves when they exist. 

Loran-C operates in the low frequency band at 100 KHz with 
a 10 KHz bandwidth. Within geographic areas, there is one 
master station and two or more secondary stations. All 
stations sequentially transmit a pulse modulated continuous 
wave in groups of eight pulses separated by one microsecond 
with the exception of the master station which transmits nine 
pulses. The ninth pulse is separated from the first eight by 
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Figure 5. Intersection Of Radio Beacon Signals Showing Error 
Spread As Range From The Source Increases. (Hobbs, 1981, 
p. 235) 



two microseconds and is used primarily for visual 
identification when viewing the pulses on an oscilloscope. 
This can be seen in Figure 7. Each group of master and 
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Figure 6. Areas Of The World Covered By Loran-C. In 
Particular, Note That The Coverage Is Almost Completely In 
The Northern Hemisphere. (Hobbs, 1981, p. 237) 



secondary stations is called a chain. Currently, there exist 
13 chains throughout the world. (USNO, 1972, p. 7-5) (Hobbs, 
1981, pp. 236-239) 

Each master and secondary station transmits according to 
a predetermined fixed order synchronized using precise atomic 
time standards. Within each receiver, the group of eight 
pulses are integrated to form one pulse of 320 microseconds. 
The time difference of arrival is then measured between the 
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Figure 7. Loran-C Pulse Sequence For A Four Station Chain. 
(Hobbs, 1981, p. 239) 

master and selected secondary station pulse from which may be 
plotted a line of position on Loran-C charts. As shown in 
Figure 8, Loran-C charts have the number coded hyperbolic 
lines overlaid on charts and to enable plotting. Chain and 
station identification are achieved by measuring the pulse 
repetition rate (PRR) and the pulse repetition interval (PRI) . 
(USNO, 1972, p. 7-5) (Hobbs, 1981, pp. 236-239) 

Accuracy of the system varies from about 700 ft near the 
baseline to about 2000 ft near the extreme range of the 
system. Due to the low frequency, baseline distances range 
from 1000 to 1500 miles so that coverage within a chain 
includes a great deal of area. Similar to the use of radio 
beacons, Loran-C charts and tables are required for the signal 
information to be meaningful. 
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Northeast U.S. Chain. (GRI 9960) (Hobbs, 1981, 242) 

D. DECCA NAVIGATOR SYSTEM 

Decca was conceived by an American and developed by the 
British Admiralty Signals Establishment during World War II 
for use in mine sweeping and the Normandy invasion forces in 
1944 (Beck, 1971, p. 65) . Today, it is operated by the Racal- 
Decca Navigator Company in only a few but very heavily 
traveled areas in the world as shown in Figure 9. It operates 
using the principles of hyperbolic radio navigation but on a 
much smaller scale than Loran-C. 

Each separate area of coverage has a master station and up 
to three slave stations arranged in a star pattern around the 
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p. 64) 

master station and identified by a color designator of green, 
red or purple. Each station within master/slave chain 
transmits an integer multiple of an unmodulated, continuous 
wave base frequency (/) in the range of 14.00 KHz to 14.33 
KHz. The master station transmits at 6/, the purple station at 
5/, the red station at 8/ and the green station at 9f. Each 
master/slave pair is phase locked along the baseline. A Decca 
receiver consists of four receivers, one for each of the 
frequencies transmitted by the master and slave stations. 
(Beck, 1971, pp. 65-67) 

Within the receiver, each master/color frequency is 
integer multiplied to get another unique multiple of the base 
frequency which is the first common multiple of either of the 
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two. For instance, the master station frequency of 6/ when 
multiplied by four and the red station frequency of 8/ when 
multiplied by three produce a new harmonic of 24/ for the 
master/red station pair. Each half wave-length of this new 
harmonic is considered a "Decca lane", and by comparing phases 
of the two signals, one can determine position within the lane 
as Figure 10 illustrates. 




Signals Along The Base Line. (Hobbs, 1981, p. 251) 

Lane ambiguities are resolved once each twenty seconds 
when all four stations transmit the fundamental frequency / 
whose half-wavelength precisely encompasses the common 
multiples (18/ for the master/green pair, 24/ for the 
master/red pair and 30/ for the master/purple pair) . By using 
phase comparison of the received signal in each of these 
"Decca zones" as shown in Figure 11, the exact lane can be 
determined. In a like manner, zone ambiguities are resolved by 
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transmitting on a frequency of 8.2/ also every 20 seconds. 
Called an orange frequency, it forms a coarse hyperbolic 
pattern within which each 360° difference cycle encompasses 




Figure 11. Phase Comparison Of Decca Master/Green-Slave Lane 
Identification Signals. (Hobbs, 1981, p. 252) 



The effective range of each chain is approximately 240 
miles from the master station. At 100 miles, the accuracy is 
approximately 30 yds by day and 100 yds by night. Again, this 
system requires the Decca charts and receiver in order to 
provide positioning data. (Bowditch, 1984, pp . 1058-1059 ) 

E. CONSOL 

Consol was developed in 1945 from a German radio beacon 
system called Sonne. It is a hyperbolic radio navigation 
system with a very short baseline such that a "degenerated" 
hyperbolic pattern is formed. The curved portions of the 
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hyperbolic pattern are not used, and so it acts as a long 
range radio beacon. It is used mostly in Europe and the 
Mediterranean. Another system very similar to Consol, called 
Cosolan, is used in the United States. Both have the advantage 
that no special equipment is required other than a low 
frequency receiver. The operator only needs to listen to the 
number of dots or dashes to know where within a particular 
sector the receiver is located. (Beck, 1971, 113-116) 

Both systems use low frequency radio waves in the range 
from 190 KHz to 350 KHz. Each station consists of three radio 
towers spaced three wavelengths apart. One of the towers 
transmits a continuous wave frequency while the other two 
transmit waves that undergo 180° phase shifts with respect to 
the CW wave during a keying cycle. All the signals are 
modulated with dots and dashes. This arrangement produces a 
series of sectors 10° to 15° wide as shown in Figure 12. The 
phase shift of the combined signal results in a variable 
number of dots and dashes being heard depending on the radial 
position within the sector. (Beck, 1971, 113-116) 

Consol range can be up to 1500 miles but not less than 25 
miles due to the compressed width between sectors close to the 
transmitters. Keeping in mind that only radial lines of 
position may be obtained from a Consol station, accuracies of 
less than a degree may be obtained at the maximum range 
improving as the receiver approaches the station. Other than 
knowledge of the station location and the dot and dash 
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7-29) 

position correlations, little is required to use Consol. 
(USNO, 1972, p. 7-5) 

F. VOR/DME AND TACAN 

VOR/DME stands for VHF Omni-directional Ranging/ Direction 
Measuring Equipment and TACAN stands for Tactical Air 
Navigation. It is, in effect, a radio beacon system for 
aircraft along heavily travelled airline routes and between 
military airfields and ships. The principle of operation of 
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TACAN is the same as for VOR/DME except that TACAN uses higher 
frequencies and is more accurate (Beck, 1971, pp. 26-28) . Some 
VOR stations are equipped to respond to TACAN DME frequencies 
and are called VORTAC stations. VOR/DME uses VHF frequencies 
between 108 MHz and 118 MHz and TACAN uses UHF frequencies 
between 960 MHz and 1,215 MHz (Logsdon, 1992, p. 36). 
Consequently, both systems are line of sight only. 

VOR employs two different signals working together in 
partnership: a narrow-beam rotating "lighthouse" transmission 
coupled with an "blinking" omni-directional pulse. As the 
lighthouse transmission sweeps past magnetic north, the 
blinking pulse is transmitted. Receiver equipment marks the 
time of receipt for the first pulse and then the time of 
receipt for the lighthouse sweep. Since the lighthouse sweep 
rotates at a fixed 30 rev/sec, it is a simple matter to 
compute the bearing from magnetic north. VOR provides only 
bearing information. (Logsdon, 1992, pp. 36-38) 

For range from the VOR station, the aircraft transmits a 
series of interrogation pulses to the DME at the VOR station 
which then transmits a reply a fixed time after receipt. 
Equipment on the aircraft computes the slant range by dividing 
the time from interrogation to the time of receipt of the DME 
response by two and multiplying by the speed of light. Figure 
13 illustrates this operation. (Logsdon, 1992, pp. 36-38) 

At the extreme range of receipt, VOR can be off by as much 
as three degrees. However, as the range to the transmitter 
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decreases, the error is reduced. Overall, VOR/DME TACAN 
positioning accuracies range from approximately 200 to 600 
feet (Logsdon, 1992, p. 36). The obvious advantage over a 
radio beacon system is lack of need for RDF equipment on board 
the aircraft. 

G. OMEGA 

The Omega navigation system was designed to provide 
worldwide coverage and to be accurate. It consists of only 
eight stations, each using a base frequency of 10.2 KHz, and 
transmitting at 10 Kw power. This makes it possible to receive 
at least three and usually four stations anywhere on earth 
with an accuracy of approximately one nautical mile. Omega's 
principle of operation is similar to Decca in that it is a 
hyperbolic radio navigation system that uses the phase 
comparison of two CW transmissions to obtain a line of 
position. It is different in that any two stations can be used 
to obtain a line of position and with special receiving 
equipment, range from the two stations for a fix. (Bowditch, 
1984, pp. 1016-1036) 

Two 10.2 KHz continuous waves transmitted exactly in phase 
but travelling in opposite directions produce a series of 
hyperbolic lanes called Omega lanes. Each lane is eight 
nautical miles wide which corresponds to the half-wavelength 
of the frequency. Phase comparison determines the receiver 
position within the lane which is expressed in terms of 
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Figure 13. VOR/DME Principles Of Operation. (Logsdon, 1992, 



P-37) 



centicycles (cec) or centilanes (cel) . These are defined to be 
0.01 of the width of the lane. Figure 14 shows two Omega lanes 
produced by the 10.2 KHz frequency. (Hobbs, 1981, pp . 259-262) 
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PHASE DIFFERENCE 

Figure 14. Two Half-Wavelength Omega Lanes Produced By The 
Phase Comparison Of The 10.2 KHz Signals. (Hobbs, 1981, p. 
260) 

Lane ambiguity is resolved by multiplexing the basic 
frequency with three others at 11.05 KHz, 11.33 KHz and 13.6 
KHz over a ten second period according to the commutation 
pattern shown in Figure 15. Within the receiver, the 
difference frequencies of 0.283 KHz, 1.133 KHz and 3.4 KHz are 
extracted which are 1/36, 1/9 and 1/3, the frequency of the 
10.2 KHz signal, respectively. The 3.4 KHz signal is used to 
establish a broader phase lane that encompasses three Omega 
lanes. Likewise, the 1.133 KHz signal phase lane encompasses 
three of the 3.4 KHz phase lanes and so on. By performing 
successive phase comparisons as demonstrated in Figure 16, the 
operator can easily determine within which Omega lane his 
receiver is located. (Hobbs, 1981, pp. 259-262) 



24 



No* way (A] 



Librru |B) 



Hawaii (C) 



N Dakota CD) 



La Reunion {£) 



A 


B 


C 


D 


! 


E 


' 1 


G 


H 




















10 2 Lf 


13.6 


U 111/3 LJ 


i2. r 


u 


12 . 1 * |_ 


J " 05 U 


12 . 1* 


U ~ 


1 2 .oi_r 


10.2 


U 136 U 


11 1/3 


U 


12.0* l 


J 1 2 . C*LT 


11.05 


U 1 2 . 0 * U 


-rr?u 


1 1 .0 


in^ir 


13.6 


u 


n 1/3 [ 


JTTTLT 


11 8 


TJ 11 - 05 u 


11.06 i_r 


U.f [J 13.1* U 


10 2 


u 


136 [ 


J 11 1/3 LI 


13. 1 4 


L! u..*U 


12.3U 


11 05 


Lh*.3 4 U 


12.3 


*u 


10.2 [ 


J ,3 - 6 Lf 


11 1/3 


UI2.3-U 



A»gtniin» IF) | 12.$^ [ 1 ? . g- | | 11.05 | | 1 2 . 9 * [_] 12.9’ |J 10.2 |J i7$ |J 1. 1/3 [J 

11 1/3 [J 1 3.0* U 1 3 .0* U ' ii.05 U 13.0* U 13.0*|J iol [j U.6 (J 



Australia IG) 
Ana 



Japan f H I 



13 6 [J 1 » V3 y 12.C* [J 1?.B 4 |J ~~n 05 [j 1 2 . jf [JTcT^ ' J 102 [J 

h 

L 



p- 0.9^1 p-l.Q-H p- 1.1-^ p— 12-^: p— 1.1— * : p-Q.9*H p—1.2 — H p-I.O—i | 



0.0 1.1 2.3 3.6 5.0 6.3 

Time Iseconds) 

4 is unique IreQu^ncy ai each station 



1A 



8 8 



100 



Figure 15. Omega Station Commutation Pattern. (Hobbs, 1981, 

p. 262) 



H. SHIP ' S INERTIAL NAVIGATION SYSTEM (SINS) 

Although it requires no external input, SINS is still an 

extremely useful navigation system. It is found primarily on 

submarines and aircraft carriers. Inertial navigation is... 

"... the process of directing the movements of a vessel 
based on sensed accelerations in known spatial directions 
by means of instruments that mechanize the Newtonian laws 
of motion, integrating such accelerations to determine 
velocity and position." (Hobbs, 1981, p. 281) 

In essence, an inertial navigation system such as SINS is a 

highly accurate dead reckoning computer. 
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Resolution On Surface Ships. (Hobbs, 1981, p. 261) 

SINS mounts two accelerometers on a platform stabilized by 
a system of three gyros in such a way that it is constantly 
coincident with a plane tangential to the earth's surface. The 
two accelerometers are continually oriented in a north-south 
and east-west direction; hence, they are sensitive only to 
horizontal north-south and east-west accelerations. 
Integrating the accelerations with respect to time gives the 
component velocities from which the ship's true velocity may 
be computed. (Hobbs, 1981, p.281-282) 

Because of the cumulative affect of a number of possible 
sources for errors such as friction in the gyro supports and 
disturbances caused by the daily rotation of the earth, SINS 
accuracy degrades over time and must be updated from external 
sources. One improvement to SINS has been the Electro-Static 
Gyro (ESG) which is used to reset or update the SINS. Simply, 
the ESG is a one centimeter diameter beryllium sphere spinning 
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at 216,000 rpm in a near perfect vacuum. This rotor is 
supported solely by an electrostatic field. The sphere is thus 
freed from the classical gyro bearing friction as well as many 
of the associated random torques that a mechanical support can 
introduce. The ESG extends the interim period between external 
updates by a factor of six. (Kayton, 1990, pp. 194-206) 

Another improvement in SINS technology is the Ring Laser 
Gyro (RLG) . It differs from the traditional concept of a gyro 
in that it contains no spinning mass. It is a closed geometric 
laser path (usually triangular such as the one in Figure 17) 
centered on an expected spin axis. Identically phased laser 
beams are continuously generated, which travel in opposite 
directions around the closed path. Any rotation about the spin 
axis causes an apparent phase difference in the two beams at 
a measurement point at one vertex of the triangle, 
proportional to the speed of rotation. When combined with the 
output of three mutually orthogonal accelerometers, position 
and velocity are easily computed. RLG's increase the precision 
over an ESG monitored SINS by an order of magnitude. (Hobbs, 
1981, p.281-282) 

Finally, in an effort to further improve accuracy and 
reduce size and weight, the latest development is the Fiber 
Optic Gyro. Its principle of operation the same as RLG's 
except the laser path can be tightly wound in a spool and 
mirrors are not required. The technology is very young and 
promising and is still being perfected. 
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Figure 17. Single Axis Ring-Laser Gyro. By Monitoring The 
Interference Pattern Caused By The Two Counter Directional 
Isophase Light Beams, Extremely Small Rotations Can Be 
Measured. (Logsdon, 1992, p. 107) 



I. SHIPBOARD DOPPLER SONAR SYSTEM 

Another navigation system that does not provide absolute 
position as much as it provides accurate inputs to a dead 
reckoning system is the Shipboard Doppler Sonar System. 
Because of its extreme accuracy in measuring very slow 
velocities, it is found on very large ships where ship's speed 
over ground is important. For example, when mooring to a pier, 
the ship could crush the pier if the contact velocity were 
greater than a knot (Hobbs, 1981, p. 285). 

It operates by measuring the doppler shift of two or more 
sonar beams in order to compute speed over ground. One such 
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system by Raytheon, called Janus, uses two fixed hull mounted 
transducers that transmit pulsed sonar signals in two beams, 
one forward and the other aft, at about a 30° angle from the 
vertical. Figure 18 illustrates a similar configuration that 
uses four transducers. The former of the two produces fore and 
aft speed accuracies of 0.1 kt, depth accuracies to the 
nearest foot, meter or fathom and distance accuracies to the 
nearest hundredth of a nautical mile. (Kayton, 1990, pp. 188- 
193) 

J. TRANSIT 

Transit Navigation Satellite System, sometimes called just 
Transit or SatNav, can trace its conceptual origins directly 
to the launch of Sputnik I. It consists of five satellites in 
polar orbits at an altitude of 580 nm with periods of 
revolution of 107 minutes as shown in Figure 19. The rest of 
the system includes ground tracking stations, a computing 
center, an injection station and receivers. Its geometry is 
configured such that every satellite comes within range of 
every position on earth at least twice a day at 12 hour 
intervals. The average time between fix opportunities 
approximately 100 minutes at the equator and a half an hour 
near the poles. The difference is due to the intrinsic nature 
of polar orbits. (Bowditch, 1984, pp. 1066-1095) 

Transit employs two frequencies in the UHF band of 150 MHz 
and 400 MHz. It operates by using an integrated doppler method 
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Raytheon Doppler Sonar System. (Hobbs, 1981, p. 285) 
where the slant range from the receiver to the satellite is 
computed by fitting the doppler phase shift curve of the 
satellite to look up table. Each satellite transmits a phase 
modulated message that repeats every two minutes containing 
the satellite time, ephemeris and other information from which 
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Figure 19. The Distinctive "Bird Cage" Polar Orbits Of 
Transit Satellite System. (Hobbs, 1981, p.275) 



the receiver calculates the receivers position in space and 
computes the longitude and latitude. Within in each pass of 
the satellite, the receiver needs to collect at least three 
successive messages for a fix. Four are preferred and seven 
are optimal. The number receivable is dictated by the height 
above the horizon achieved by the satellite. (Bowditch, 1984, 
pp. 1066-1095) 
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The accuracy of Transit is affected by ships dead 
reckoning accuracy and speed of advance. For a stationary 
ship, it is common to achieve accuracies of 3 5 meters made 
possible by the shorter wavelength of the UHF signals 
(Bowditch, 1984, pp. 1067). It is significant to note that 
while this system is world wide and has improved positioning 
accuracies of better than an order of magnitude over systems 
discussed so far, it is comparable to Omega in its coverage 
and worse in its availability to provide positioning 
information. 

K. GLOBAL POSITIONING SYSTEM (GPS) 

With the June 1993 launch of the twenty-fourth satellite, 
NAVSTAR's Global Positioning System (GPS) became fully 
operational in late 1993 . GPS provides real time three 
dimensional positioning and timing information anywhere in the 
world. It consists of 24 uniformly spaced satellites in six 
orbital planes at an altitude of 20,180 Km as depicted in 
Figure 20. Like Transit, it uses two UHF frequencies at 1,575 
MHz and 1227 MHz simultaneously upon which constellation 
ephemeris data, atmospheric propagation correction data and 
satellite clock error are phase modulated. (Wells, 1987, pp. 
4 . 00-4 . 11) 

The principle of operation for GPS is to calculate the 
range from four different satellites using the broadcast 
message, and to compute the longitude, latitude and altitude 
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Least Twelve In View At All Times Anywhere In The World. 
(Hobbs, 1981, p. 288) 

of the simultaneous solution of the resulting ranges. Because 
the system is operated by the Department of Defense for 
military use, the lower of the two frequencies is encrypted 
and altered to deny any potential enemy real time positioning 
accuracy for medium and long range targeting. Still, 
accuracies using just the higher frequency are approximately 
100 meters (Logsdon, 1992, p.64). Using methods to be 
explained in the next chapter, this can be improved by more 
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than an order of magnitude and in non-real time applications 
such as surveying, mind boggling accuracies of a millimeter 
spherical error probable (SEP) can be achieved using phase 
comparison techniques (Logsdon, 1992, p.64). 

L. SUMMARY OF CONTEMPORARY ELECTRONIC NAVIGATION SYSTEMS 

It is easy to see why, after surveying the various 
electronic navigation systems available, GPS invites the most 
enthusiasm and excitement from all users. It is a quantum 
improvement over anything else available. There are other 
navigation systems not discussed here such as Microwave 
Landing Systems (MLS) , the French Argos system and the Russian 
GPS imitation GLONASS. All suffer from one or more limitations 
in range, coverage, availability or accuracy. What is 
important in this review is an appreciation of the various 
contemporary electronic methods and their method of operation 
including their particular limitations. In the next chapter, 
the GPS principle of operation is more closely examined as 
well as the various methods for its use. 
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GLOBAL POSITIONING SYSTEM 



GPS operational and technical characteristics must be 
understood in detail in order to follow the reasoning for GPS 
receiver selection and the concept of employment in 
interfacing them with the CPU. In the following sections, a 
comprehensive review of the User Segment of GPS is provided. 
It covers the salient points of signals, coding, message 
content, message structure, calculations, and various biases. 
The chapter ends by using this information in a description 
and discussion of Differential GPS. 

GPS is functionally divided into three segments: Control 
segment, Space segment and User segment. The Control segment 
exists to monitor, update and operate the whole system. The 
Space segment consists of the 24 satellites. These two 
critically important segments, however, are of no concern to 
the user and little is required to be known of them. The User 
segment is the most important, for understanding its operation 
will greatly influence the user's choice of receiver and 
subsequent method of employment. 

A. USER SEGMENT RECEIVER HARDWARE DESCRIPTION 

Currently, receivers cost as little as a few hundred 
dollars to as much as $50,000 or more. Size and weight also 
vary greatly from receivers as small as a package of 



35 



cigarettes weighing approximately the same to larger than a 
car battery and weighing only slightly less. Shape and type of 
antenna can also impact the effectiveness of the receiver. 
Figure 21 illustrates the most common designs of antennae 
found on the market. 

One type of antenna called a microstrip is relatively thin 
(approximately 1/10" to 1/2" thick) and can be designed to 
receive one or both of the frequencies that GPS satellites 
transmit. This design is ideal for use in this project where 
low drag is desired on an aircraft. However, signal reception 
from satellites low on the horizon usually experience a much 
lower signal-to-noise ratio due to the lower gain of the 
antenna at those low angles caused by its thin design. 

Other factors that distinguish receivers from one another 
are the choice of position solution algorithm and interface 
communication protocols. All of the above characteristics are 
determined and limited by the state of the art of the 
technology and by the purpose for which it was designed. 

Receivers may be loosely grouped into four broad 
categories based on capability (Logsdon, 1992, p. 66): 

1. Number of channels and sequencing rate. 

2. Access to selective availability signals. 

3. Use of available performance enhancement techniques. 

4. Computer processing capabilities. 

The majority of user receiver sets sold are continuous 
tracking or parallel receivers. They typically have from four 
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Figure 21. Four Common Types Of GPS Antennae. (Wells, 1987, 
p. 7.02) 



to twelve channels. There are also single channel receivers 
capable of providing relatively crude positioning. To obtain 
data from multiple satellites, they employ one of three 
different sequencing speeds: slow, fast or multiplexed (very 
fast) as shown in Figure 22. (Logsdon, 1992, p. 66-67) 
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Figure 22. Multiplexing And Sequencing Channels For GPS 
Receivers. (Wells, 1987, p. 7.13) 



Selective Availability (SA) is the primary means by which 
the United States Department of Defense denies the potential 
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terrorist or enemy accurate use of GPS for targeting 
navigation for missiles. GPS messages are encoded using two 
codes. SA alters and encrypts the message content of the more 
accurate Precision (P) code carried on both frequencies. The 
less accurate Course and Acquisition (C/A) code message is 
called the Standard Positioning Service (SPS) and is available 
to anyone on the higher frequency only. The P-code information 
is called the Precise Positioning Service (PPS) . Only those 
receivers with a decoding lip (tightly controlled by the 
DOD) can receive PPS information. (Logsdon, 1992, p. 67) 

SPS receivers can still achieve excellent accuracies by 
using performance enhancement techniques. One such technique 
capable of centimeter accuracy uses interferometry. This 
process needs very expensive antennae and receivers, and the 
data has to be post-processed. Surveyors are the most common 
users with this accuracy requirement and do not need the real 
time positioning information so denial of PPS is not a 
significant handicap. Another capability that enhances SPS 
accuracy in near real-time uses a process called Differential 
GPS (DGPS) . Simply, DGPS employs two receivers: one stationary 
at a known location and one mobile. The stationary receiver 
computes the difference between its known location and its GPS 
derived location. It then transmits those differences to the 
mobile receiver which uses them to correct its own GPS derived 
position. Since the Selective Availability errors are the same 
over baselines as large as 600 nm, accuracies approaching and 
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exceeding those of PPS receivers can be achieved. DGPS is 
explained more completely later. (Logsdon, 1992, p. 67) 

The last category of receivers is distinguished by 
computer processing capability. Mostly esoteric, this is where 
the manufacturer takes the GPS message and calculates the 
output message accounting for all the nuances of the solution 
algorithm. These include such items as whether or not to use 
just four satellites or more when they are available, whether 
or not to account for ionospheric and tropospheric delays 
(sometimes modeling these delays can introduce more errors 
than they eliminate) , when to shift to other satellites that 
will reduce the geometrical error, etc. In addition, the 
receiver can be enhanced through the use of a Kalman Filter to 
help smooth the position and time data as well as reduce 
receiver processor time. (Logsdon, 1992, p. 68) 

B. GPS SIGNALS AND CODING 

GPS satellites continuously transmit two frequencies: LI 
at 1575.42 Mhz and L2 at 1227.60 Mhz. The respective 
wavelengths for LI and L2 are approximately 20 cm and 25 cm. 
Each of these is modulated with the GPS message at 50 bps. 
Since each satellite uses the same frequencies to send its own 
unique message, the receiver must be able to distinguish one 
signal from the other. Furthermore, the transmitting power of 
each satellite so is low that by the time the signal reaches 
the receiver, its signal-to-noise ratio is less than the 
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ambient noise level. GPS overcomes these obstacles by using 
spread spectrum and code division-multiple access techniques. 
As mentioned earlier, the messages are encoded using two codes 
called C/A code and P-code. The chipping rate of the C/A code 
and P-code is 1.023 Mbps with a virtual wavelength of 290 m 
and 10.23 Mbps with a virtual wavelength of 29 m, 
respectively. Each satellite uses its own unique code with 
which the receiver may auto-correlate a duplicate and thus 
phase lock on to it. The combination of the transmission, 
coding and message chipping rates provides for a very robust 
signal with which it is very difficult to interfere. (Wells, 
1987, pp. 4.5 - 4.10) 

The C/A code uses a 1023 bit Gold code with a period of 
one millisecond. The Gold code was chosen because of its 
length and because of the property that Gold codes are nearly 
orthogonal. It is also very easy to replicate and with its 
short period, easy to achieve correlation with itself. A GPS 
signal can only be received if the code assignment is known 
ahead of time. The Gold code assigned to each satellite is 
published in the Interface Control Document ICD-GPS-200 
published by Rockwell Space Systems Division. (Wells, 1987, p. 
6.07) 

The P-code uses a much longer code with 
235,469,592,765,000 bits and has a period of approximately 267 
days. It is divided into 32 seven day segments, and a segment 
is assigned to each satellite. The DOD further encrypts the P- 
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code to produce the Y-code to deny its use to unauthorized PPS 
users. The encryption and length of the P-code make it 
extremely difficult to correlate to without keying information 
and decryption hardware. (Wells, 1987, p. 6.07) 

The GPS message is phase modulated onto its respective 
codes which in turn are phase modulated onto their respective 
carrier frequencies. As mentioned earlier, the LI frequency is 
modulated with both the C/A code and the P-code. This is 
accomplished by phase modulating its quadrature component with 
the P-code. This feature primarily allows the military to 
employ PPS redundancy in certain applications. The table in 
Figure 23 compares the important points of the code and data 
features. (Wells, 1987, pp. 6.3 - 6.8) 

C. GPS MESSAGE CONTENT 

A GPS receiver determines its position by computing the 
position of four or more satellites within view and then 
measuring its range from each. By solving four simultaneous 
equations, the receiver calculates its coordinates in an Earth 
Centered Earth Fixed ( ECEF) cartesian coordinate system. The 
system chosen is defined to be consistent with the World 
Geodetic Survey of 1984 (WGS-84) which places the origin at 
the intersection of the Z axis (through the earth's center of 
mass and parallel to the north pole of the Conventional 
Terrestrial Pole (CTP) ) and the X-axis (the intersection of 
the WGS-84 Reference Meridian Plane parallel to the zero 
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C/A Code Note: C/A code aod P code transmitted on carrier in phase quadrature 



Parameter 


C/A Code 


P Code 


Data 


CHIPPING RATE 


1.023 X 10* 
bils/sec 


10.23 X 10 fc 
bits /sec 


50 bits/sec 


SPATIAL LENGTH 
(PER BIT) 


290m 
(960 ft) 


29m 
(96 ft) 


6000km 
(3235 n. m.) 


REPETITION 

INTERVAL 


0.001 sec 


7 days 


Not applicable 


CODE TYPE 


Gold code 


267 day pseudo 
random code 


Not applicable 


TOTAL NO. OF 
CODES 


37 unique 
Gold codes 


37 seven-day sections 


Not applicable 


SPECIAL 

PROPERTIES 


Easy to 
acquire 


Slightly more accurate. 
Resistent to jamming 
and spoofing. 

Rejection of multipath. 


Provides handover from 
C/A to P code, ephcmcris 
data and clock correction. 



Figure 23. GPS Code Comparison. (Applied Research 
Laboratories, 1993, p. 2.2) 



meridian plane, and the plane of the CTP's equator). The Y- 
axis completes the ECEF right hand orthogonal coordinate 
system measured in the equatorial plane 90° to the east as 
shown in Figure 24 (Applied Research Laboratories, 1993, p. 
2.1). Having defined the reference frame for GPS position, it 
is necessary to understand a few more concepts in describing 
position within an ECEF system. 

The GPS receiver translates the ECEF coordinates into 
geodetic coordinates for the user. The equations to translate 
coordinates from ECEF to geodetic coordinates and vice-versa 
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As Defined By The WGS 84. (BIH Stands For Bureau 
International De L'Heure) (Applied Research Laboratories, 
1993, p. 2.1) 



are provided in Appendix A. By using the reference ellipsoid 
model of the world defined by WGS 84, GPS is capable of 
relating pin point positions between continents. Before the 
World Geodetic Surveys of 1972 and 1984, inter-continent 
surveying was virtually impossible. Cartographers could only 
establish local datums to which all surveyed points were 
referenced. These datums were quite accurate for their 
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respective geographical regions at best. One of the 
accomplishments of WGS 84 was in relating all the different 
local datums to the ECEF coordinate system. By modelling the 
datum as a point of tangency of an ellipsoid as shown in 
Figure 25, the ellipsoid origin could easily be referenced to 
the origin of the ECEF coordinate system. With the advent of 
WGS 84 and GPS, the importance of the local datums has been 
diminishing. Most GPS receivers contain a library of the major 
local ellipsoid references to the WGS ECEF coordinate system 
so charts and maps that use different local datums can be 
employed with GPS. 

Once the reference frame is established and understood, 
satellite position can be described using the six classical 
Keplerian orbital elements. They are briefly described in 
Figure 26 but they are not sufficient for the most accurate 
description of the satellite orbit. This is because the 
predicted position of a satellite is the result of least- 
squares curve fitting which is not valid for an entire orbit. 
Thus, it is necessary to include several other parameters to 
more accurately model satellite positions in space and time. 
All the parameters are listed in Appendix B. All told, there 
are 17 specific parameters included in the GPS message. 

D. GPS CALCULATIONS 

The equations required to determine the ECEF coordinates 
of the satellites from the ephemeris are presented in Appendix 
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Figure 25. Two Dimensional Illustration Of A Theoretical 
Ellipsoid Tangent To The Geoid At An Arbitrary Datum P. 
(Leick, 1990, p. 189) 



B. If the receivers and satellites could be perfectly 
synchronized with GPS time, three parameters would be unknown: 
X, Y, Z. It would then be possible to determine position using 
only three satellites. However, all receivers use inexpensive, 
lightweight and comparatively inaccurate quartz clocks and 
each has a bias from GPS time for which there must be an 
accounting. Therefore, for three dimensional position, there 
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THE SIX KEPLERIAN ORBITAL ELEMENTS 
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Figure 26 . Six Classical Keplarian Orbital Parameters. 
(Logsdon, 1992, p. 143) 

are four unknowns and four simultaneous equations as shown in 
Equation (2) : 



p* 1 = fix 1 - xf 1 + (y 1 - y K ) 2 + (z 1 - z K ) 2 + c ac k 

p k 2 = fix r ~ r lcf T ~^ {y 2 - r K ) 2 + (z 2 - z K ) 2 + c At K 

p K 3 = - xj 2 + (y 3 - yj 2 + (z 3 - zj 2 + cit K 

P 4 = ^/U 4 - XJ 2 + ( Y -4 - iV) 2 + (Z 4 - Zj 2 + c At, 



where P/ is the pseudorange for the Kth epoch in time, # is 
the satellite number, c is the speed of light and At, is the 
receiver clock bias from GPS time. 
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Parenthetically, it should be noted that although GPS time 
is synchronized within one microsecond of UTC, it does not 
recognize UTC leap seconds. The last simultaneous epoch of GPS 
and UTC time was in January 1980. Consequently, GPS time 
differs from UTC by the integer number of leap seconds since 
then . 

E. GPS MESSAGE STRUCTURE 

The entire GPS message format is illustrated in Figures 27 
through 29. It is defined by ICD-GPS-200 written by Rockwell 
International who was the prime contractor for GPS. The basic 
message unit is called a Frame. It is subdivided into five 
subframes which are in turn subdivided into ten words each. A 
word has 30 bits. For positioning, a receiver requires only 
subframes one through three. Those contain the necessary 
ephemeris and clock corrections. However, the designers saw a 
need to provide the user information about all the satellites 
in the entire GPS constellation. They used subframes four and 
five to do this. Since this "Almanac" of extra information 
exceeded the size of those two subframes alone, it is 
transmitted using subframes four and five of 25 successive 
message frames. Each pair of subframes four and five of the 25 
successive frames is referred to as a page. The entire block 
of 25 pages and subframes one through three is considered a 
Master Frame. Figures 27 and 28 illustrate this best. At a 
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data rate of 50 bps, a frame needs 30 seconds and a master 
frame requires 12.5 minutes for complete transmission. 



BASIC MESSAGE UNIT IS ONE FRAME (1500 BETS LONG) 

◄ 30 sec 



1 FRAME = 5 SUBFRAMES 



1 



1 SUBFRAME = 10 WORDS I 2 I 3 1 4 I 5 I 6 I 7 | 8 j 9 1 10 



6 sec 




= 30 BITS 



0.6 sec 



ZEE m 1 Hj |m itttt 



■ 0.02 sec 



ONE MASTER FRAME INCLUDES 



ALL 25 PAGES OF SUBFRAMES 4 & 5 
= 37,500 BITS TAKING 12.5 MINUTES 



Subframes 4 and 5 
have 25 PAGES 



J 



Figure 27. GPS Message Format Breakdown. (Rockwell, 1987, p. 
7.13) 



F. GPS BIAS 

There are numerous factors that affect GPS performance 
accuracy. Selective Availability is primary among these. While 
the exact nature of Selective Availability is classified, it 
is sufficient to know that it denies single receiver position 
accuracy of better than 100 meters 95% of the time. In 
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addition, the Ionosphere and Troposphere affect the time of 
receipt by slowing down the signal. If the receiver could 
receive two frequencies, these atmospheric errors could be 
determined and eliminated. For the receivers that are capable 
of only single frequency reception (most receivers) , 
compensation for these errors is achieved by using 
mathematical models which adequately reduce their affects for 
most positioning purposes. 



The GPS Data Message Content 



1 


2 


3 1 




/ \ n 


H 4 


5 


/ Orbit parameters 







Flags (L2 code & data; week #; 
satellite accuracy and health) 

Age of data 

Satellite clock correction coefficients 



Almanac for satellites 1-24 (pages 1-24) 
Health of satellites 1-24 (page 25) 



Almanac for satellites 25-32 (pages 2, 3, 4, 5, 7, 8, 9, 10) 
Ionospheric model, and UTC data (page 18) 

Antispoof flag - 32 satellites (page 25) 

Satellite configuration - 32 satellites (page 25) 

Health of satellites 25-32 (page 25) 



Figure 28. GPS Message Content Organization Within The 
Message Format. (Applied Research Laboratories, 1993, p. 
2 . 2 ) 



Signal reflections, commonly referred to as multipath, are 
another source of position errors. Multipath interferes with 
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the fundamental concept of measuring the time from 
transmission to the time of receipt in order to determine 
range from the satellites. Fortunately, with adequate 
filtering and by paying careful attention in choosing the 
antenna, almost all the effects of multipath may be 
eliminated. 



The GPS Data Message 

50 Bits Per Second Data Rate 




* Format of Frame Twenty-Five Changes. 

•* ) 2.5 Minutes Before The Entire Message Repeats. 

Figure 29 GPS Message Organization. (Applied Research 
Laboratories, 1993, p. 2.2) 

Finally, satellite geometry relative to the receiver plays 
a role second only to Selective Availability in determining 
the precision achieved. The measure of "goodness" of satellite 
geometry is called Geometric Dilution Of Precision or GDOP. 
Optimal GDOP (the lower the GDOP number the better) is when 
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three of the satellites are on the horizon equally spaced 120° 
apart and the fourth is directly over head as shown in Figure 
30. Since the constellation is designed such that there are 
usually 12 satellites in view at all times, a low GDOP of two 
or three can regularly be achieved. 




Positioning 

accuracy 



= DOP • a 

t 

Geometry 
(Dilution of 
Precision) 



Measurement 
accuracy 





POOR GDOP 
satellites bunched 
together 



GOOD GDOP 
(ideal case) 
one satellite overhead 
3 on horizon, 

120° apart in azimuth 



Figure 30. Geometric Relation Between GPS Satellites Affects 
The Position Accuracy. (Wells, 1987, p. 4.22) 



G. DIFFERENTIAL GPS 

SPS does not guarantee accuracies less than 300 m. It is 
possible to achieve accuracies two orders of magnitude better 
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through a process called Differential GPS but it requires at 
least two receivers. The concept is simple and generates the 
greatest amount of promise in commercial applications of GPS. 
By placing one of the two receivers in a pre-surveyed location 
where its geodetic coordinates are known as precisely as 
possible, it can compare its true position to the position GPS 
provides and calculate the differences between latitude, 
longitude and altitude. If the resulting error signs are 
changed, the errors become corrections that may added to the 
second receiver's GPS position for a greatly improved position 
accuracy. Over baselines of up to 600 nm, the errors are 
approximately the same. 

In practice, there are a few restrictions in employing 
DGPS. The most important is that the position results from 
both receivers must be derived using the same satellites or 
else the improvements in position are diminished by the 
introduction of new errors. The other restriction is that the 
corrections should be applied at the earliest point in the 
algorithm used to determine position coordinates in order to 
minimize the propagation of errors. The reasons for these 
restrictions will become clear after a brief description of 
the algorithm used to compute DGPS corrections. 

The stationary GPS receiver is receiver is required to 
send corrections to the mobile receiver. It would be simple to 
subtract the true longitude, latitude and altitude from the 
geodetic position obtained from the GPS satellites. However, 
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the errors introduced from SA, satellite geometry and 
atmospheric interference affect the psuedoranges from each 
satellite in varying degrees. One satellite psuedorange will 
always have a greater error than the next. Since the GPS 
position is calculated using a least squares reduction of the 
four simultaneous equations listed in equation two, all the 
individual psuedorange errors will be "blended" thereby 
permitting the greatest error to act on all the psuedoranges 
and the final position solution. 

In order to avoid this problem, the stationary receiver 
computes the expected position of and range from each 
satellite (based on the time of travel of the individual 
signal) and compares it to the position and range that the GPS 
message asserts in its ephemeris. The difference between the 
two becomes a psuedorange correction. At least four of these 
corrections are transmitted to the mobile receiver and applied 
to the corresponding psuedoranges obtained by that receiver. 
Then the new psuedoranges are used in equation two to 
calculate the longitude, latitude and altitude of the mobile 
receiver. Since the individual corrections are applied before 
equation two, the final position is more accurate. 

In practice, GPS makes all of its calculations using the 
Cartesian Coordinate ECEF system so the stationary receiver 
must translate its true position from geodetic coordinates to 
ECEF coordinates. In a similar manner, the mobile receiver 
must translate its corrected position from ECEF coordinates to 
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geodetic coordinates. To summarize, the following steps are 
involved in improving the GPS position using DGPS: 

1. Convert the true geodetic position of the stationary GPS 
receiver to ECEF coordinates. 

2. Compute the satellite positions from the ephemeris. 

3. Compute the expected satellite positions by multiplying 
the time from transmission to time of receipt by the speed 
of light. 

4. Compare the received pseudoranges to the expected 
pseudoranges and calculate the pseudorange corrections for 
each satellite. 

5. Transmit the pseudorange corrections to the mobile GPS 
receiver and add to the received pseudoranges. 

6. Compute the mobile receiver's position using the 
corrected pseudoranges. 

7. Convert the mobile receiver's position from ECEF 
coordinates to geodetic coordinates. 

H. SUMMARY 

The User segment of GPS draws together the results of 
several decades of study and experimentation. Hardware design 
and miniaturization are pushing the state of the art. Signal 
construction and frequency selection optimize its immunity to 
noise and interference. Orbital kinematics referenced to an 
ECEF coordinate system simplify and standardize the 
information contained in the GPS message structure. Even so, 
natural and man-made biases corrupt the integrity and accuracy 
of GPS. Through DGPS users can overcome the degradations 
caused by these biases and achieve real time accuracies 
unheard of heretofore. Understanding these principles upon 
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which the User segment is founded will enable the user to 
better select a receiver that suits his purpose and integrate 
it into whatever application he so desires. 
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IV. GPS RECEIVER SELECTION AND INTERFACE WITH CPU 



A. MISSION REQUIREMENTS FOR GPS 

Before any serious search for a GPS receiver can be made, 
it is necessary to define the extent to which it will be used. 
These criteria come from the mission statement. This project's 
goals are to integrate an Inertial Measurement Unit (IMU) with 
GPS and other non-INS sensors through a smoothing and control 
filter to achieve autonomous flight of an Unmanned Aerial 
Vehicle using off the shelf hardware and technology. As a 
result, several requirements are defined for the GPS 
receivers . 

To begin with, the position output should be frequent 
enough to provide regular corrections to the IMU. Also, the 
GPS receiver must actually update its position solution as 
often as possible to ensure accuracy in the information of its 
outputs. The receiver that is to be in the aircraft should be 
small and lightweight. It should be easily interfaced to the 
CPU, and it should provide the necessary information to 
conduct phase measurements of the incoming GPS signal. (This 
last feature should provide the positioning accuracy to auto- 
land the aircraft.) Finally, it must be readily available on 
the open market. Intrinsic in the mission statement is that it 
must be affordable given the project budget. With these 
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guidelines, an informed selection can be made in choosing the 
best GPS receivers for the job. 

B. GPS RECEIVER SELECTION 

Having reviewed the principles of GPS operation and 
variations in receiver design features, an extensive survey of 
216 different GPS receivers from 54 companies was conducted. 
In the January 1993 issue of GPS World magazine, many of the 
critical facets of each receiver were listed for comparison 
(Chan, 1993, pp. 52-65). In order to stay within budget, only 
those receivers less than $5,000.00 were considered. (The 
highest priced receiver was $55,000.00 with many between 
$10,000.00 to $20,000.00.) Among the remaining units, those 
that were capable of tracking six or more satellites 
simultaneously were selected for comparison. This left the 
list at seventeen receivers. 

Of the remaining seventeen receivers, all were capable of 
the same approximate accuracy and information output rate at 
once per second. The one notable exception was Navstar's XR5 
with an output rate at 1/10 second, but ultimately, it was not 
chosen because, at $1000.00, it cost three times as much as 
the Motorola PVT-6. Since all the receiver units performed at 
virtually the same level, the deciding factors came down to 
cost, weight, size and power consumption. Table I shows the 
top eight receivers with these features listed. 
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Table I. COMPARISON OF TOP EIGHT CANDIDATES FOR THIS 
PROJECT'S GPS RECEIVERS. (CHAN, 1993, PP. 52-65) 



GPS 

Maker 


RCVR 

NAME 


Max # 
chan 


Weight 

j 


Size (in 3 ) 


Power 

(Watts) 


Price 


Furuno 


GB-92/ 

GN-72 


8 


2.8 


8.74 


2 


$750 


Magellan 


GPS 

Brain 


5 


4.8 


6.53 


H . 3 


$445 


Magnavox 


Turbo Eng 


8 * 


3.2 


11.44 


<1.5 


$3750 


Motorola 


PVT-6 


8 . 


4.5 


4.64 


1.3 


$369 


Navstar 


XR5 

Kernel 


12 


5 


8.28 


2.5 


$995 


Novatel 


Perform 

OEM 


10 


6.2 


15.44 


5 


$3195 


Sony 


IPS- 

2010 


8 


8 


15.65 


2.61 


$1000 


Trimble 


SV6 


6 


2.9 


2.85 


1.85 


$995 



The Trimble price quoted includes a starter evaluation 
kit. The Motorola receiver, when all options were installed, 
was actually more expensive. By the criteria used, Trimble had 
the better unit would have been chosen over Motorola but they 
were very slow in answering inquiries for characteristics and 
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pricing information. When the decision was made, they were not 
considered. However, two factors make the choice of Motorola's 
PVT-6 the best selection. One, they also make a portable 
differential GPS base station called a LGT-1000. It is 
conceivable, in the future, this project will choose to 
experiment with the specialized equipment. And two, Trimble 
has since stated they do not intend to pursue the small, 
single receiver user market. Instead, they will specialize in 
industrial scale receivers for use at airports and in 
surveying. 

C. GPS RECEIVER DESCRIPTION 
l. Hardware 

The Motorola PVT-6 is an original equipment for 
manufacturer (OEM) module. Although evaluation software is 
provided in a starter kit, it is not useful when the CPU is 
multi-tasking with several other processes. PVT-6 stands for 
Position, Velocity and Time with six channels. Figure 31 shows 
the dimensions of the unit without the antenna. The 
specifications are listed in Appendix C. The PVT-6 is capable 
of tracking six satellites simultaneously. It uses only the 
C/A coded message from the LI carrier. The code tracking is 
carrier aided. It can be powered by a 12 Vdc unregulated or 5 
Vdc regulated source. (Motorola, 1992, p. 3-2) 

The PVT-6 is supplied with a low profile antenna of a 
microstrip design called a strip patch antenna. The antenna 
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Figure 31. Physical Dimensions Of The PVT-6 . (Motorola , 1992, 

P- 3-3) 



gain pattern is provided in Appendix D. The LI frequency is 
collected by the antenna, passed through a narrow band filter, 
amplified and sent to the receiver unit via an RG-58 coaxial 
cable. It provides the antenna 5 Vdc @ 25 mA from the receiver 
and returns the amplified LI signal from the antenna. It 
should be noted the power loss along the cable should not 
exceed 6 dB at a frequency of 1575.42 MHz which equates to a 
maximum length of six meters for RG-58 coaxial cable. RG-58 is 
the recommended cable because it allows a sufficient length of 
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cable while achieving lighter weight due to the small gauge 
without the corresponding losses due to increased resistance 
that accompany even smaller gauge coaxial cables at that 
frequency. Antenna placement should position the patch plane 
to be as level as possible and have an unobstructed view of 
the horizon. This ensures a direct line of sight to all the 
visible satellites. Table II and Figure 32 provide the power 
and data connector pin assignments. (Motorola, 1992, p. D-l) 

Table II SERIAL PORT AND POWER PIN ASSIGNMENTS. (MOTOROLA, 



1992, p. 3-2) 



Pin it 


Signal Name 


Description 


1 


+12V/+5V BATT 


(Optional) +5 Vdc regulated or +12 Vdc 
unregulated for running Real Time Clock and 
retention of satellite ephemeris information 
stored in Static RAM memory. 


2 


+5V MAIN 


(Optional) +5 Vdc regulated for power 
requirements of the entire GPS receiver(+12 Vdc 
signal not used). 


3 


+12V/+5V RTN 


Power supply (+5V or +12V) return. 


4 


V 1 

PP 


Flash memory (EPROM) programming voltage. 


5 


+12V MAIN 


+12Vdc unregulated for power requirements of the 
entire GPS receiver. 


6 


1 PPS 2 


(Option A) 1 pulse per second output. 


t 


1PPS RTN 


(Option A) 1 pulse per second return. 


8 


RS232 TXD 


Serial RS232 data output. 


6 


RS232 RXD 


Serial RS232 data input. 


10 


RS232 RTN 


Signal return for RS232 signals. 



'Not used (used only by Motorola for updating of software) . 

: 1 Hz pulse train with rising edge coincident with UTC/GPS 1 second 
tick. 



As Figure 33 shows, the LI analog signal enters the 
receiver and is down converted to an intermediate frequency 
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ANTENNA DATA/POWER 

CONNECTOR CONNECTOR 




Figure 32. Side View Of The PVT-6 Showing The Pin Numbering 
Order. (Motorola, 1992, p. 3-2) 

(IF) . The IF is passed to the 6-channel code and carrier 
correlator section where it is digitized and split for input 
into six separate channels for correlation, filtering, carrier 
tracking, code tracking and signal detection. The output is 
synchronously routed to the MPU (PVT-6 CPU) where the 
information is processed and formatted. Also within the MPU is 
the RS-232 interface for full-duplex, asynchronous data 
communication with the host equipment. (Motorola, 1992, p. 2- 
2) 

The PVT-6 is configured to operate at 9600 b~ id with 
eight bit words, one start and stop bit and no parity. Its 
communications port is considered Data Communications 
Equipment with respect to the RS-232 serial communications 
standard. However, close examination of the available pins 
for serial communications reveals the presence of only the TX 
(pin 8) , Rx (pin 9) and Gnd (pin 10) pins. This dictates there 



63 




Figure 33. Functional Block Diagram Of The PVT-6. (Motorola, 
1992, p. 2-1) 



will be no modem line configuration (DTR and DSR) and at least 
no hardware flow control (CTS and RTS) . In fact, there is no 
software flow control (XON and XOFF) either implying the Data 
Terminal Equipment (DTE- host equipment) must monitor the port 
continuously or miss valuable information. Likewise, the PVT-6 
must do the same. In practice, it receives all serial input 
into a two kilobyte buffer from which it reads 750 bytes per 
second. The buffer is large enough that it never overflows. 
(Motorola, 1992, pp. 3-2 - 4-3) 
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Finally, the PVT-6 is provided with static random 
access memory (SRAM) for retention of the satellite ephemeris 
data. To prevent the loss of this information, an external 
power source such as a battery may be hooked up to the 
+12V/+5V BATT connection (pin 1) when the GPS unit is powered 
off. In addition, nonvolatile EEPROM is used for storage of 
custom operating parameters, almanac information and other 
information such as the receiver serial number and option 
list. (Motorola, 1992, p. 2-2) 

2. Motorola Binary Message Format 

The PVT-6 receiver is capable of providing its 
information in any one of three formats: Loran, NMEA-0183 and 
Motorola Proprietary Binary (Motorola, 1992, pp. 4-2 - 4-6). 
Both Loran and NMEA (National Marine Electronics Association) 
would severely restrict this project's flexibility and 
accuracy, and so are not used. Motorola Binary format is 
logically organized and provides tremendous flexibility to the 
user to apply GPS to what ever purpose. 

Motorola's Binary format has 47 different messages it 
sends and receives ranging in length from 7 to 294 bytes. Each 
message begins with "@@" (double ASCII 0x40) which is followed 
by a two letter, case sensitive message ID. Then, 0 to 287 
bytes of data follow, a byte size check sum and an ASCII 
carriage return-line feed (OxOD OxOA) . The check sum is the 
exclusive-or of all the preceding bytes in the message. The 
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data that follows the message ID can be thought of as scaled 
floating point data stored in integer format. This avoids the 
floating point storage convention as well as saves space (only 
one byte is needed instead of four or more) . The user must be 
very careful to observe the units of measurement provided in 
the Motorola GPS Technical Reference Manual. (Motorola, 1992, 
p. 4-2) 

For every command to the GPS receiver, there is a 
response usually using the same message ID. The receiver has 
no way to determine if the information it has received from 
the CPU is correct other than by verifying the several 
intrinsic properties of the message itself: 

• The message begins with 

• The message ends with CR-LF. 

• The length of the message matches that expected for the 
message ID. 

• The checksum is valid. 

• The message data is within the expected maximum and 
minimum ranges. 

If the command input to the PVT-6 changing a parameter is 
invalid for any of the above reasons, it responds as expected, 
but the data values will not reflect the ordered change. It is 
prudent, therefore, to compare every input message with the 
subsequent PVT-6 response. (Motorola, 1992, pp. 4-2 - 4-3) 

Input commands may be of the type that change a 
particular configuration parameter of the receiver. Any 
command with a message ID that begins with a upper case "A" 



66 



falls into this category. Input commands may also be of the 
type that enable or disable the output of data or status 
messages. They are the commands that the CPU will use for 
measuring position, velocity, time, pseudorange and satellite 
ephemeris data. These output status messages include 12 that 
are the work horses of the PVT-6. Table III shows only nine 
because three of the commands are the reciprocals used by the 
CPU to talk to the PVT-6. The nine shown are listed in order 
of precedence for transmission from the GPS receiver to the 
CPU. (Motorola, 1992, p. 4-3) 

Table III DATA MESSAGE OUTPUT RATES SHOWN IN ORDER OF 
PRECEDENCE OF TRANSMISSION .ROM THE PVT-6 TO THE CPU. 
(MOTOROLA, 1992, P. 4-4) 



|{!|; OUTPUT- MESSAGE TYPE jj&jy 


1$ . j 




Position/Channel Status 


At selected update rate 


When requested 


Satellite Range Data Output 


At selected update rate 


When requested 


Pseudorange Correction Output 


At selected update rate 


When requested 


Ephemeris Data Output 


When Eph data changes 


When requested 


Visible Satellite Status 


When Vis data changes 


When requested 


DOP Table Status 


When DOP data changes 


When requested 


Almanac Status 


When Aim status changes 


When requested 


Almanac Data Output 


When Aim data changes 


When requested 


Leap Second Pending 




When Requested 



3. PVT-6 Operation 

The PVT-6 has some limitations that, if understood, 
may be overcome. Data is output once per second from the 
receiver. This information is calculated from the ephemeris 
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and time measurements of the previous second time tic called 
the measurement epoch. In order to provide the most accurate 
position given the latency of data, the position, velocity and 
time represent the best estimate of the data advanced one 
second. Figure 34 graphically illustrates this process. The 
first byte of information for a particular measurement epoch 
is transmitted anywhere from zero to 50 ms after the 
subsequent measurement epoch. Only the time for the Satellite 
Range Data message and the Pseudorange Correction message is 
not propagated forward. Careful consideration addressed to 
this latency issue in the design of any Kalman Filter or 
smoothing filter that integrates GPS position with the IMU and 
other sensors will improve performance. (Motorola, 1992, pp. 
4-13 - 4-17) 

To assist in synchronization of the GPS data with 
other sensor data, a one pulse per second (1PPS) signal is 
output using pins six and seven of the serial port. The pulse 
is approximately 200 milliseconds long with its leading edge 
corresponding precisely to the measurement epoch. The PVT-6 
has the option of additionally offsetting the 1PPS up to 
0.999999999 seconds. If there is a requirement for the 1PPS to 
be precisely synchronized with the UTC or GPS measurement 
epoch, the user may also offset the receiver measurement epoch 
as well as account for antenna cable delay. For this project, 
these latter conveniences will not be required. This is shown 
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Figure 34. Position/Channel/Status Message Data Output 
Latency. (Motorola, 1992, p.4-15) 



in the timing diagram in Figure 35. (Motorola, 1992, pp. 4-16 
- 4-17) 



Finally, the PVT-6 suffers from periodic outages like 
most other receivers regardless of quality. These outages may 
be caused by temporary loss of satellite track do to obscuring 
by a building, tree or wing tip if the aircraft were in a 
banking turn. Sudden accelerations may also cause loss of 
track. In this event, the time for the receiver to reacquire 
the track and output data is important. As shown in Table IV, 
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Figure 35 Timing Diagram For The Epoch, 1PPS And Serial 
Output Data. (Motorola, 1992, p. 4-17) 



this time may take as long as 15 seconds. If the aircraft were 
attempting to land using GPS alone, this could be fatal. 
(Motorola, 1992, p. 4-7) 

Table IV NOMINAL REACQUISITION TIMES FOR THE PVT-6. 

(MOTOROLA, 1992, P. 4-7) 



OBSCURED 


IfltB^GCRJfSFl^^ 




15 sec 


< 12 sec 


30 sec 


<13 sec 


45 sec 


< 1 4 sec 


60 sec 


<15 sec 
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D. INTERFACING THE GPS RECEIVERS 

The two GPS receivers to be used for DGPS will be linked 
to the ground CPU and the UAV CPU through RS-232 serial ports. 
The big picture of the project showing where GPS fits in is 
shown in Figure 36. The ground control station will be 
interfaced through a 33 MHz Intel 80486DX personal computer. 
It will be the primary data flow conduit for the uplink and 
downlink as well as provide formatting and mass storage. The 
additional link in the ground station shown in Figure 36 is 
for redundancy in the event the main link fails and emergency 
control of the UAV is required. 




Figure 36 Archytas Project Control/Communications Diagram 
Showing Where GPS Fits Into The Big Picture. 
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The UAV CPU uses off the shelf hardware that economizes on 



size, weight and price. This is done by using six computer 
cards and a passive backplane. The CPU is a 33 MHz Intel 
80486DX with 256 Kb double cache memory. Non-volatile memory 
for the operating system and programs is provided by a 2.88 Mb 
RAM/ROM Disk card instead of a relatively heavy hard drive. In 
addition, servo control and sensing are managed using an 
Analog-to-Digital ( A/D) /Digital-to-Analog (D/A) card and a ten 
channel Timer/Counter card. Serial I/O is handled through the 
two standard serial communication ports on the CPU card (C0M1 
and COM2) and an intelligent 8 port RS-232 serial I/O card. 
The GPS receiver communicates to the CPU through port 1 on 
this card (not to be confused with C0M1 on the CPU card) . In 
order save space, all the cards are connected to a common bus 
using a passive backplane. A simplified drawing of the final 
arrangement is shown in Figure 37. 

1. PCL-744 8-Port Intelligent RS-232 Interface Card 

The UAV CPU needs multiple serial ports for other data 
inputs from the link modem, the IMU and the non-INS sensors. 
These multiple serial inputs are provided through Advantech's 
PCL-744 8-Port Intelligent RS-232 Interface Card. 

The PCL-744 is an intelligent card in that it has its 
own processor and 64 Kb RAM buffer. It requires only one 
interrupt number for eight ports although each may configured 
independently. This configuration includes determining the 
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size of the transmit and receive ring buffers and what type of 
flow control is used. These facts enable the UAV CPU to 
perform multi-tasking by allowing it to access the ports and 
process data without losing any information from the other 
ports arriving simultaneously. The PCL-744 takes care of that 
by multiplexing between its eight ports and placing the 
incoming data in buffers that can be accessed by the UAV CPU. 

The PCL-744 comes with a library of 40 user functions 
written in several high level languages including Microsoft C 
and Turbo C by Borland. This project uses the Turbo C routines 
for all serial communications. These functions allow the user 
to use the card in anyway desired for serial I/O. It is also 
necessary to ensure the TSR drivers are loaded prior to PCL- 
744 operation. They are called STD-DRV.EXE and 744-DRV.EXE and 
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should reside in the AUTOEXEC.BAT file to ensure their 



loading. 

The description of the software drivers for the GPS 
receivers that follows is the same in functional design for 
both the ground and UAV units, but there are slight 
differences in coding. This is due to the fact that the UAV 
receiver must communicate through the PCL-744, and the ground 
receiver communicates directly with the CPU through a standard 
DOS COM port. The port interface for the ground CPU uses a 
standard C library interrupt routine, and the port interface 
for the UAV that communicates through the PCL-744 uses the PC- 
ComLIB serial communication library routines provided with the 
card . 

2. GPS Software Driver 

In communicating with the PVT-6, all transmissions are 
without any flow control. It is necessary, therefore, to 
provide a storage buffer where the contents of the 
transmission may be stored for holding until the CPU is ready 
to handle the data. The information sent to the PVT-6 is put 
into a two kilobyte buffer in the receiver itself. All message 
handling is done by the receiver including parsing, 
verification and interpretation, and is transparent to the 
user. The user's only responsibility is to ensure that the 
message sent to the receiver is properly formatted. However, 
when receiving messages from the PVT-6, parsing, verification 



74 



and interpretation must be done by the user. This allows 
flexibility for the user and is one of the minor reasons for 
selecting a Motorola GPS receiver. 

This section describes the software driver required to 
accept information from the GPS receivers. The software in 
Appendix E is written in modular format to allow flexibility 
and ease of handling. It is structured so it can be expanded 
and built upon. The following paragraphs describe the flow 
chart shown in Figure 38 which illustrate the concept of the 
GPS software driver in Appendix E and the structure of the 
software driver. 

The PCL-744 requires the eight ports to be initially 
configured using the vendors software. This configuration is 
stored in the PCL-744 driver and the card is thus reconfigured 
upon startup every time power is applied. This does not make 
the card ready for operation, though. Each port must be 
individually opened and set to the communication parameters 
required by its user. For instance, the PVT-6 uses 9600 BAUD 
with an eight bit word, no parity and one start/stop bit. In 
addition, RS-232 modem control (i.e., DTR and RTS), hardware 
control (i.e., DTS and RTS) and software flow control (i.e., 
XON/XOFF) are set to off. This is all accomplished using the 
function "Open_ports" . At the conclusion of using the ports 
they should be formally closed and this is accomplished by 
using the function "Close_ports" . 
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Figure 38 Functional Diagram Of The GPS Software Drivers For 
Receiving Information From The PVT-6 . 



The main functions that handle incoming information 
are called "Master_gps" for the ground CPU and "Slave_gps" for 
the UAV CPU. The "Master_gps" function first gets the buffer 
length of the message waiting in the Rx buffer for port number 
one on the PCL-744. It then creates a buffer of exactly that 
size and extracts the information from the Rx buffer and puts 
it into the temporary buffer just created. Next, it creates 
the input buffer which will contain the temporary buffer 
appended to the excess buffer contents. The excess buffer has 
the partial message left over from the previous call to 
"Master_gps" . With the input buffer created and the contents 
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of the excess and temporary buffers copied into it, the excess 
and temporary buffer memory is freed. 

The pointer to the input buffer is then passed to the 
function "Parse_inbuf f " where the first message is parsed from 
the input buffer and returned. If a complete message is not 
returned, the function places the partial message in the 
static excess buffer, frees the memory for the input buffer 
and exits the function "Master_gps" . If a complete message is 
returned, it is passed to the function "Msg_verif ication" 
where it is checked for the correct length, check sum, ASCII 
prefix and hexadecimal suffix xOD xOA (ASCII carriage 
return and line feed) . When it has been verified, the message 
is placed in a buffer to be passed out to the calling routine 
for transmission to the UAV GPS if it is a pseudorange 
correction. Otherwise, it is written to disk for post- 
processing . 

In the UAV GPS driver, the "Slave_gps" function 
operates exactly the same way except after a message has been 
verified, it is put into a structure and scaled according to 
the format in the Motorola protocol for subsequent use by the 
Kalman filter. One of the more complicated structures for the 
ephemeris message has been written as an example by Motorola 
and incorporated in the header file GPSCONV.H. All message 
output from the UAV GPS is passed out to the calling routine 
to be down linked to the ground in binary format (to reduce 
downlink time) for post-processing. 
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Two other functions are mentioned but not yet 
completely written which play an important role. They are the 
"Initialize_master_gps" and "Initialize_slave_gps" functions. 
They will pre-conf igure the receivers for such things as true 
geodetic coordinates of the ground antenna, time, date, output 
rates, desired files to be collected etc. They should be made 
interactive to allow the user to change and set any 
parameters. The code framework for these functions is in the 
file GPSPORTS . H in Appendix E. 

Finally, the source file GPSAN.C with GPSCRNCH.H was 
used to post-process data collected for this thesis. It was 
written to take the Position/Channel/Status message from an 
ASCII file, convert it to binary and then place the message in 
a structure where the scaled and formatted data could be 
analyzed. Then GPSCRNCH.H "crunched" the numbers and computed 
the maximum, mean and standard deviation for the various time 
delays. 
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V. VERIFICATION OF OPERATIONAL CAPABILITIES 



GPS is a complicated new tool for navigators. Each 
receiver manufacturer employs unique and innovative techniques 
to improve its accuracy and flexibility. Motorola provides 
specifications and operational capabilities for its receiver 
and users have published various accuracies achieved using 
DGPS. It is prudent to verify those features important to this 
project's success. Those two capabilities are examined below. 
They are: 

1. Verification of the manufacturer's reacquisition time 
specifications . 

2. Determination of the effects of pseudorange correction 
delay in getting from the stationary receiver to the mobile 
receiver . 

In the following sections, these concerns are investigated. 
For each case, the experimental setup is described and the 
results are presented. 

A. REACQUISITION TIME VERIFICATION 

This project intends to use DGPS with one of the two 
antennae on a UAV. During some flight maneuvers, it is 
expected that the GPS antenna will be masked from one or more 
of the satellites it is using to calculate its position. It is 
important for the designers of the Kalman Filter to understand 
the normal reacquisition times given such a circumstance. In 
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Table IV in Chapter IV, Motorola presents its results for the 
reacquisition time of the PVT-6 receiver given the antenna was 
masked for 15, 30, 45 and 60 seconds. The method they used to 
collect this data is not well documented, and it is not 
indicated whether or not the receiver was in three dimensional 
fix mode or two dimensional fix mode. The difference could be 
significant because the former requires four satellites for a 
fix and the latter only three. 

The reacquisition time was investigated using this 
project's PVT-6 receiver. Table V shows the results next to 
the Motorola claims. Data was collected using the software 
supplied by Motorola with an evaluation kit. The antenna was 
completely masked 360° around above and below using four to 
six layers of tin foil. The mask was left in place for the 
specified amount of time and then removed for at least 30 
seconds. Data files were imported into a spread sheet where 
the time marked data was analyzed. 

The results do not exactly agree with those that Motorola 
publishes, but they are close. Significantly, data obtained 
during daylight hours were within one standard deviation of 
Motorola's claimed reacquisition times whereas data collected 
after daylight hours were not. No effort was made to confirm 
this suggested phenomenon nor was any effort made to correlate 
it to ionospheric thickness or tropospheric conditions such as 
rain or cloudiness. However, until the IMU is integrated with 
GPS using a Kalman Filter and the effects of a two to three 
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second delay in reacquisition time can be determined, 

investigation of this discovery is left for further study. 

Table V COMPARISON OF MOTOROLA AND EMPIRICAL REACQUISITION 
TIMES. 



Time Obscured 

(sec) 


Motorola 
Reacquisition 
Time (sec) 


Empirical 
Reacquisition 
Time (sec) 


Standard 

Deviation 

(sec) 


15 


<12 


13 . 75 * 


2.9 


45 


<13 


17 . 08 ** 


2.0 


45 


<14 


16 . 67 ** 


3 . 2 


60 


<15 


13 . 67 * 


2.1 



* These observations were taken from 1230 to 1500. 
** These observations were taken from 2200 to 2400. 



B. DGPS ACCURACY AND PSUEDORANGE CORRECTION LATENCY 

The accuracy that DGPS can achieve is well documented. The 
PVT-6 is designed and built to permit DGPS. In fact, the 
Motorola Binary format includes master/slave DGPS commands and 
automatically incorporates the data in its calculations. It is 
worthwhile to validate the expected results especially in 
light of the fact the communications software to be used is 
not written by the manufacturers and may adversely affect the 
accuracy . 

In designing the test for DGPS, the length of the baseline 
between the two receivers was an issue that was researched. 
The results indicated the baseline range between the two 
receivers did not significantly affect DGPS accuracy up to 
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about 50 km. Beyond 50 km, out to about 1000 km, the effects 
slowly increased the error (Kremer, 1990, pp. 307-312). This 
project's maximum baseline is not expected to exceed 15 km. 
For all intents and purposes, a zero meter base line is as 
good as a 1000 meter baseline and certainly is easier to 
test. Therefore, this test simply mounted the two antennae 
side by side at a pre-surveyed location (See Appendix F) . 
Data was collected by sending pseudorange corrections from the 
master antenna to the slave antenna using the software driver 
in Appendix E. The PVT-6 puts out a position message at a rate 
of once per second. These positions were differenced using the 
ECEF equations in Appendix A, and the mean, maximum and 
standard deviation were calculated. 

Table VI shows the DGPS accuracy with various delays in 
transmitting the correction to the slave receiver and Figure 
39 graphically illustrates the same results. The delay was 
artificially introduced through the function call to 
"Pr_delay" in the header file GPSDELAY.H. The delays were in 
integer amounts of seconds as shown below. 

Interestingly, Motorola technicians quote adverse effects 
of pseudorange corrections received after 90 seconds. The 
results above suggest an immediate degrading effect of any 
amount of delay and that the error grows with delay time. It 
is hypothesized that the decrease in error at the 15 second 
delay point is a spurious result caused by any number of 
factors. It is left for further study to confirm this. 
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Table VI OBSERVATION OF PSEUDORANGE CORRECTION INPUT DELAY 
ON THE EFFECTIVE ACCURACY OF DIFFERENTIAL GPS. 



Delay 


Average 


Maximum 


Standard 


(sec) 


Offset 


Offset 


Deviation 




(m) 


(m) 


(m) 


0 


1.17 


3.76 


1.00 


1 


2 . 58 


5.16 


0.87 


• 


2.75 


3 . 1 ? 6 


1 . 00 


3 


2 . 72 


3 . 86 


0.96 


4 


0.96 


3.28 


0.92 


5 


2.41 


5.53 


0.97 


6 


1.85 


3.85 


1.39 


7 


2.45 


8.03 


2 . 02 


8 


4 . 07 


8 . 83 


2 . 14 


9 


4 . 51 


10.80 


2 . 19 


10 


4 . 94 


13.85 


3 .36 


15 


1.96 


5.91 


1.48 


30 


6.09 


21.21 


4 . 59 


45 


21.40 


42.16 


10.90 
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Showing The Effects Of Time On The Accuracy Of DGPS 
Corrections . 
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VI. CONCLUSION 



This thesis reviews the major contemporary electronic 
navigation systems available and demonstrates that GPS offers 
the best accuracies and world wide coverage on a realtime 
basis. The user segment of GPS is explained, and its 
capabilities and deficiencies discussed. Differential GPS is 
presented and shown to be able to achieve enhanced accuracies 
sufficient to significantly contribute to this project. 

An extensive review of available GPS receivers was 
conducted and Motorola's PVT-6 was found to combine all the 
benefits of GPS in a small, lightweight, low power and 
inexpensive package. Motorola's proprietary binary interface 
protocol allows the user great flexibility with little 
sacrifice of convenience. In addition, the standardized RS-232 
serial communications port enables easy access to a variety of 
different systems. All things considered, it is the best 
choice over all other GPS receivers for this project. 

The PVT-6 's reacquisition times compare favorably with the 
manufacturer's claims. When two of these receivers were used 
in a DGPS mode with a zero baseline, the results easily 
confirmed nominal results of one to three meters desired 
accuracy. The most encouraging result is the relative immunity 
of DGPS accuracy to short delays (one to seven seconds) of the 
pseudorange correction inputs. 
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The success of this project's ultimate goal: to perform an 
autonomous landing, will depend on improving processing speeds 
and accuracies. A few areas where further investigation into 
GPS could contribute to this end are listed below. 

• Conduct dynamic tests to confirm expected nominal 
accuracies using DGPS. 

• Conduct additional data collection to determine if the 
accuracy of DGPS or reacquisition time is affected by the 
time of day or by the weather. 

• Collect more data to determine if psuedorange correction 
latency steadily degrades or has anomalies (wherein the 
distance error actually improves with time) as Figure 39 
suggests at the 15 second mark. 

• Determine the effects of the aircraft's RF environment on 
the GPS performance. 

• Investigate whether programming in assembly or some other 
language might improve processing speeds. 

• Determine if the PVT-6's capability to compute two 
dimensional position by constraining altitude (relying 
solely on barometric altitude sensor data for altitude 
inputs) might improve robustness of GPS satellite 
tracking . 

• Write a serial port driver that is independent of the PCL- 
744 multiple serial port I/O card in order to provide 
flexibility in using the PVT-6 and possibly reduce the 
total hardware weight through the elimination of the PCL- 
744 itself. 

• Write a mission initialization interface routine to setup 
up the PVT-6 receivers. 

• Write a post-processing routine to interpret and present 
the recorded results in a desired format. 

• Write a real time graphical interface to plot the 
coordinates on a map. Also include either a separate 
window or color coded routine to be able to see current 
altitude . 
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GPS has revolutionized electronic navigation. The airline 
industry and other nations decry the U.S. DOD's introduction 
of Selective Availability. They are actively lobbying the U.S. 
government to force the DOD to turn it off and leave it off. 
It is an uphill battle but not a hopeless one. In time, as 
civilian applications become common place in our society, the 
DOD may acquiesce and turn it off. If that happens, GPS 
accuracies achieved in this project will undoubtedly improve 
and may warrant reconsideration in the way they are used in 
the Kalman Filter that integrates data from all the sources. 
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APPENDIX A: COORDINATE TRANSFORMATIONS 



A. GEODETIC COORDINATE TRANSFORMATION TO ECEF COORDINATES 

Given the geodetic coordinates longitude (X) , latitude (0) 
and altitude (h) , the Earth Centered Earth Fixed coordinates 
( X,Y,Z ) are (Leick, 1990, p. 184-185): 

X - ( N + h ) cos<}) cosX 
Y = { N + h ) cos<}) sin! 

Z = [N(l-e 2 ) + h ] sin<J) 

where h is the height above the ellipsoid and 

e 2 = 2 f - f 2 

~ _ a - b 
a 

a = semi -major axis 
b = semi -minor axis 

N = ____ a 

/~1 - e 2 sin0 



B. ECEF COORDINATE TRANSFORMATION TO GEODETIC COORDINATES 

Given the ECEF coordinates ( X,Y,Z ), the Geodetic 

coordinates (X,0,h) are (Leick, 1990, p. 184-185) : 
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0 = arctan 



1 - 



X = arctan _ 
X 



h = 

COS0 



Z 

e2 \i~ X 2 + X 2 



- N 
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APPENDIX B: EPHEMERIS ALGORITHM 



The following parameters are provided continuously for 
each satellite of the GPS constellation. Their exact placement 
and description within the ephemeris message is defined by the 
Interface Control Document, ICD-GPS-200 as published by 
Rockwell International Space Systems Division. Any serious 
manipulation of the ephemeris algorithm using GPS data 
requires careful consideration of this document. There are 
several important, yet subtle, facts about the time tagging 
and the corresponding validity of the Issue of Data messages. 



M 0 


Mean Anomaly at Reference Time 


An 


Mean Motion Difference from Computed Value 


e 


Eccentricity 


(A) 1/2 


Square Root of the Semi-Major Axis 


fl 0 


Longitude of Ascending Node of Orbit Plane at 
Weekly Epoch 


io 


Inclination Angle at Reference Time 


w 


Argument of Perigee 


OMEGADOT 


Rate of Right Ascension 


IDOT 


Rate of Inclination Angle 


Cue 


Amplitude of the Cosine Harmonic Correction 
Term to the Argument of Latitude 


C U s 


Amplitude of the Sine Harmonic Correction Term 
to the Argument of Latitude 
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C rc Amplitude of the Cosine Harmonic Correction 

Term to the Orbit Radius 

C rs Amplitude of the Sine Harmonic Correction Term 

to the Orbit Radius 

C ic Amplitude of the Cosine Harmonic Correction 

Term to the Angle of Inclination 

C is Amplitude of the Cosine Harmonic Correction 

Term to the Angle of Inclination 

t w Reference Time Ephemeris 

IODE Issue of Data (Ephemeris) 

The following tables list the number of bits, the location 
within words three through ten (omitting the parity bits) , the 
scale factor of the LSB, the range and the units for each of 
subframes two (Table VII) and three (Table VIII) . Parameters 
indicated with an asterisk are to be two's complimented with 
the sign bit in the MSB. Unless otherwise listed in the value 
range column, effective range is the maximum range attainable 
with the indicated bit allocation and scale factor. 

The following elements of the ephemeris algorithm are 
listed in order of expected solution. All the required 
information is given or may be found using the above from the 
GPS ephemeris message. 

WGS 84 value of the earth's 
H = 3.986005 x 10 14 meters universal gravitational 

sec 2 parameter 
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Table VII BIT ALLOCATION AND SCALING FACTORS FOR WORDS THREE 
THROUGH TEN OF SUBFRAME TWO. (ROCKWELL, 1987, PP. 60-81) 



Parameter 


No. of 
bits 


Location 
of bits 


LSB scale 
factor 


Value 

range 


Units 


IODE 


8 


1-8 






see ICD— 
GPS-200 


c rs 


16* 


9-24 


2 s 




meters 


An 


16* 


25-40 


2 43 




semi- 

circles/ 

sec 


Mo 


32* 


41-72 


2-31 




semi- 

circles 


^UC 


16* 


73-88 


2 29 




radians 


e 


32 


89-120 


2-33 


0.03 


dimension 

less 


C u , 


16* 


121-136 


2- 29 




radians 


(A) 1/2 


32 


137-168 


2-1 9 




meters 172 


toe 


16 


169-184 


2 4 


604,784 


seconds 


other 


8 


185-192 






see ICD- 
GPS-200 



*Parameters indicated are signed integers. All others are unsigned. 



, WGS 84 value of the earth's 

0 = 7.2921151467 x 10' 5 rotation rate 

sec 



A = (v/3) 2 



Semi-major axis 



n 



o 




rad 

sec 



Computed mean motion 



t k = t 



'oe 



Time from ephemeris reference 
epoch 



In the above equation, t is GPS system time at time of 
transmission (i.e., GPS time corrected for transit time 
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Table VIII BIT ALLOCATION AND SCALING FACTORS FOR WORDS 
THREE THROUGH TEN OF SUBFRAME THREE. (ROCKWELL, 1987, PP 
60-81) 



Parameter 


No. of 
bits 


Bit 

location 


LSB scale 
factor 


Value 

range 


Units 

l 


C, c 


16* 


1-16 


2' 29 




radians 


n 0 


32* 


17-48 


2 31 




semi- 

circles 


C 1S 


16* 


49-64 


2-29 




radians 


i-o 


32* 


65-96 


2-31 




semi- 

circles 


C rc 


16* 


97-112 


2 5 




meters 


0 ) 


32* 


113-144 


2‘ 31 




semi- 

circles 


OMEGADOT 


24* 


145-168 


2- 13 




semi- 

circles/ 

sec 


IODE 


8 


169-176 






see ICD- 
GPS-200 


I DOT 


14* 


177-190 


2 ^ 




semi- 

circles/ 

sec 


other 

! 


2 

1 . J 


191-192 

L _ __ J 


1 




see ICD- 
GPS-200 , 



*Parameters indicated are signed integers. All others are unsigned. 



range/speed of light. . Furthermore, t k shall be the actual 
total time difference between the time t and the epoch time 
t^., and must account for beginning or end of week crossovers. 
That is, if t k is greater than 302,400 seconds, subtract 
604,800 seconds from t k . If t k is less than 302,400 seconds, 
add 604,800 seconds to t k . 



n = n 0 + An 



Corrected mean motion 



= M 0 + nt k 



Mean anomaly 
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The next equation solves for the eccentric anomaly (E k ) by 
iteration. M k and e are known from above. 



Let 



Then 



E, = 



M 



k 



1 - e 



M, = E, - e sin E 



A, 



M k - M, 

1 - e cos E 



If M 2 ~ M k , then E 2 = E k 
otherwise repeat and 
compare again. 

(Tattersf ield, 1984, p. 36) 



E 2 - E, + A, 



Having found E k , the rest of the solution is simply solving 
equations . 



v 



k 



= arctan 



sfl ~ e 2 sinE k 
cosE k - e 



True anomaly 



$ = v k + co 



Argument of latitude 



Second harmonic perturbations: 
6u k =C us sin20 k + C uc cos2$ k 



6 r k =C rc cos 2<E> k + C rs sin2<E> k 
Si k =c ic cos20> k + C is sin2«D k 

u k = $ k + 5r k 

r k = A (1 - e cosE k ) + 5r k 
i k = i 0 + Si k + ( IDOT) t k 



Argument of lat correction 
Radius correction 
Correction to inclination 

Corrected argument of latitude 
Corrected radius 
Corrected inclination 
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Positions in orbital plane 



Xk = r k cosu k 
y k = r k sinu k 



^0 + ^ ) ^k ^e^oe 



Corrected longitude of 
ascending node 



x k = x k cosQ k 
x k = x k sinQ k 
z k = Yk sini k 



Yk cos i k sinQ k 
y k cosi k cosQ k h 



Earth-Centered, Earth-Fixed 
(ECEF) coordinates 
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APPENDIX C: PVT -6 SPECIFICATIONS 
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APPENDIX D: PVT- 6 ANTENNA GAIN 
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APPENDIX E: SOFTWARE DRIVER SOURCE CODE 



A. GPSAN.C 



r GPSAN.C 

Reads ascii position data "@@Ba n from a file, converts it 
it to binary and places the scaled and formatted messages 
in an array of structures. 

Input is the name of the ascii file without the extension. 

Output is a file with the analyzed data from GPSCRNCH.H. 

*/ 

^include <string.h> 

^include < stdlib. h > 

^include <stdio.h> 

^include "c:\borlandc\twite\gpstruct.h" 

^include M c:\borlandc\twite\gpsconv.h M 
^include "c:\borlandc\twite\gpscrnch.h" 

^include " c:\borlandc\twite\gpsfun.h" 

void main(void) 

{ 

^define COL 130 
^define NEWCOL 66 

char name(40], path[30] = "c:\\borlandc\\etwork\\"; 
int i,j,k,m,n,s, master, slave, ROW; 

ONEBYTE first, second, third, inbufftCOL]; 

ONEBYTE * *mbuff, **sbuff, **mbuffbin, **sbuffbin; 
struct T_POS_CHAN_STATUS *marst, *sarst, pos_chan; 

FILE *fpname; 

printf("\nEnter the input filename (without the extension(??.ext): "); 

scanf("%s",name); 

strcat(path, name); 

strcatlpath, ”.txt”); 

printf("%s",path); 

fpname = fopen(path,"r"); 

printf("\nEnter the number of rows (individual messages): ”); 

scanf("%d",&ROW); 

printf("\n"); 
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k = 0; 
m = 0; 
s = 0; 

/* Create mbuff and sbuff. 

*/ 

mbuff = calloc< ROW, sizeofl ONEBYTE * )); 
sbuff = callocf ROW, sizeof( ONEBYTE * )); 
for ( i = 0; i<ROW; i + + ) 

{ 

mbuffli] = calloc( COL, sizeof( ONEBYTE )); 
sbuff [i] = calloc( COL, sizeof{ ONEBYTE )); 

}; 



while (1 ) 

{ 

I* Read in file to ascii arrays mbuff and sbuff. 

*/ 

if ( k > = 70 ) 
break; 

first = fgetc(fpname); 
if ( first = = ) 

{ 

second = fgetc(fpname); 
if (second = = ) 

{ 

third = fgetc(fpname); 

I* Store inbuff array in either mbuff or sbuff. 

*/ 

if (third = = 

{ 

fread( sbuff [s], sizeof(ONEBYTE), COL, fpname ); 
k = 0; 
s + + ; 

} 

else 

{ 

mbufftmHOl = third; 
for (i = 1 ; i<COL; i + + ) 

mbuff[m](i] = fgetc(fpname); 
k = 0; 
m + + ; 

}; 

}; 

}; 

k + + ; 

}; I* end while * l 
fclose(fpname); 

master = m; 
slave = s; 
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r Create mbuffbin and sbuffbin. 

*/ 

mbuffbin = callocf master, sizeoff ONEBYTE * )); 
for ( i = 0; i< master; i + +) 

mbuffbin[i] = callocf NEWCOL, sizeoff ONEBYTE )); 
sbuffbin = callocf slave, sizeoff ONEBYTE * )); 
for ( i = 0; i< slave; i + + ) 

sbuffbinfi] = callocf NEWCOL, sizeoff ONEBYTE )); 

I* Convert hex buffer to binary buffer. 

*/ 

for (i = 0; i < master; i+ + ) 

{ 

mbuffbin[i][0] = mbuff[i][0]; 
mbuffbin[i][1] = mbuff[i][1 ]; 
k = 2; 

for ( j = 2; j < COL; j+ =2) 

{ 

mbuffbin[i][k] = (Hex(mbuff[i][j]) < < 4) | Hex(mbuff[i][j+ 1 ]); 
k + = 1; 

}; 



for (i = 0; i < slave; i + + ) 

{ 

sbuffbin[i][0] = sbuff [i] [0]; 
sbuffbin[i][1 ] = sbuff [i][ 1 ]; 
k = 2; 

for ( j = 2; j < COL; j + = 2) 

{ 

sbuffbin [i][k] = (Hexfsbuf f [ i] [j]) < < 4) j Hexfsbuff [i] [j + 1 ]); 
k + = 1; 

}; 

}; 



I* Free mbuff and sbuff memory. 

*/ 

for (i = 0; i<ROW; i + + ) 

{ 

free(mbuff[i]); 
f reefsbuff [i]) ; 

}; 

free(mbuff); 

free(sbuff); 

/* Temporarily store mbuffbin in a file "mtemp.bin" and free memory. 
*/ 

fpname = fopenf "mtemp.bin", "wb" ); 
for (i = 0; i < master; i + + ) 

fwrite(mbuffbin[i], sizeof(ONEBYTE), NEWCOL, fpname); 
for (i = 0; i< master; i + + ) 
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free(mbuffbin[i]); 

freefmbuffbin); 

fclose(fpname); 

I* Temporarily store sbuffbin in a file "stemp.bin" and free memory. 
*/ 

fpname = fopenf "stemp.bin", "wb" ); 
for (i = 0; i <slave; i + + ) 

fwrite(sbuffbin[i], sizeof(ONEBYTE), NEWCOL, fpname); 
for (i = 0; i< slave; i + + ) 
f reefsbuff bin[i]) ; 
free(sbuffbin); 
fclose(fpname); 

r Create an array of structures. 

*/ 

marst = callocf master, sizeoff struct T_POS_CHAN_STATUS )); 
sarst = callocf slave , sizeoff struct T_POS_CHAN_STATUS )); 

fpname = fopenf "mtemp.bin", ”rb"); 
for ( i = 0; i < master; i + + ) 

{ 

freadfinbuff, sizeof(ONEBYTE), NEWCOL, fpname); 
Convert_Pos_Chan_Status( inbuff ); 
marstfij = pos^chan; 

}; 

fclosef fpname ); 

removef” mtemp.bin"); 

fpname = fopenf "stemp.bin", ”rb"); 

for ( i = 0; i < slave; i + + ) 

{ 

freadfinbuff, sizeof(ONEBYTE), NEWCOL, fpname); 
Convert_Pos_Chan_Status( inbuff ); 
sarstfi] = pos_chan; 

}; 

fclosef fpname ); 
removef "stemp.bin"); 

Ave_Max_Stdev(marst, sarst, master); 

}; r end main */ 
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B. 



GPSCONV.H 



r GPSCONV.H 

Convert_Ephemeris() 

Convert_Pos_Chan_Status() 

Convert_xDOP() 

Convert_Sat_Range_Data() 

Convert_Pos_Chan_Status_Ext() 

*/ 

#ifndef _GPSCONV_H 
^define _GPSCONV_H 

^include <math.h> 

#ifndef _GPSTRUCT_H 

^include "c:\borlandc\twite\gpstruct.h" 

#endif 

/************* + + ♦ + *** + ******* + * + ************* + + ***********♦ + ■*■ + ** 



Convert_Ephemeris() 

This procedure converts the raw binary ephemeris data 
into a representation usable for post-processing 
applications. The equations coded in this procedure 
are taken from the ICD-GPS-200 (the GPS NAVSTAR ICD). 

Input is an array with the binary ephemeris message. 

Although it is unsealed, it is in the array 
format of the ICD-GPS-200 format. 

Output is a structure with the scaled and formatted 
ephemeris message. 

*/ 

void Convert_Ephemeris( ONEBYTE hex_ephemeris[][3] ) 

{ 

/* Sign extends a neg hex byte into its TWOBYTE (integer) equivalent 
*/ 

^define MACRO_SET_UPPER(word_to_set) \ 

if ( word to set & 0x0080 ) word_to_set | = OxFFOO ; 

I* This MACRO takes two hex bytes and merges them into a single 
TWOBYTE (integer) value. 

*/ 

^define MACRO_MAKE_EPH_WORD(e_wd, wl, bl, w2, b2) \ 

e_wd = ( ( ( UNSIGNED_TWOBYTE ) hex_ephemeris[w1 Hbl ] << 8 ) j \ 
( UNSIGNED_TWOBYTE ) hex_ephemeris[w2][b2] ); 
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UNSIGNED_TWOBYTE wordl, word2; /* Temporary storage for words built. */ 
struct T_FLOAT_EPHEMERIS *eph ; 

eph = malloc(sizeof( struct T_FLOAT_EPHEMERIS) ); 

/* Time of ephemeris 
*/ 

MACR0_MAKE_EPH_W0RD(word1 , 15, 0, 15, 1 ); 

eph- > toe = (double) ( <(UNSIGNED_FOURBYTE) wordl) * 16); 

/* Mean anomaly at reference time ( mO ). 

*/ 

MACRO_MAKE_EPH_WORD(word1 , 9, 2, 10, 0 ); 
MACRO_MAKE_EPH_WORD(word2, 10, 1, 10, 2 ); 
eph->m0 = ( ( ( FOURBYTE ) wordl < < 16 ) j word2 ) * PITWOM31; 

/* Mean motion difference from computed value ( delta_n ). 

*/ 

MACRO_MAKE_EPH_WORD(word1 , 9, 0, 9, 1 }; 
eph->delta_n = ( TWOBYTE ) wordl * PITWOM43; 

/* Eccentricity ( e ). 

*/ 

M ACRO_MAKE_EPH_WORD(word 1 , 11, 2, 12, 0 ); 
MACRO_MAKE_EPH_WORD(word2, 12, 1, 12, 2 ); 

eph->e = ( ( ( UNSIGNED_F0URBYTE ) wordl < < 16) j word2 ) * TWOM33; 

/* Square root of the semi-major axis ( sqrt_a ). 

*/ 

MACRO_MAKE_EPH_WORD(word1 , 13, 2, 14, 0 ); 
MACRO_MAKE_EPH_WORD(word2, 14, 1, 14, 2 ); 

eph->sqrt_a = ( ( ( UNSIGNED_FOURBYTE ) wordl << 16 } j word2 ) * TWOM19; 

/* Longitude of ascending mode of orbit plane at 
weekly epoch ( omega_0 ). 

*/ 

MACRO_MAKE_EPH_WORD(word1, 16, 2, 17, 0 ); 
MACRO_MAKE_EPH_WORD(word2, 17, 1, 17, 2 ); 

eph->omega_0 = ( ( ( FOURBYTE ) wordl < < 16 ) j word2 ) * PITW0M31; 

/* Inclination angle at reference time ( iO ). 

*/ 

MACRO_MAKE_EPH_WORD(word1 , 18, 2, 19, 0 >; 
MACRO_MAKE_EPH_WORD(word2, 19, 1, 19, 2 ); 
eph->i0 = ( ( ( FOURBYTE ) wordl << 16 ) | word2 ) * PITW0M31; 

I* Argument of perigee ( w ). 

*/ 

M ACRO_M AKE_EPH_WORD(word 1 , 20, 2, 21, 0 ); 
MACRO_MAKE_EPH_WORD(word2, 21, 1, 21, 2 ); 
eph-> w = ( ( ( FOURBYTE ) wordl < < 16 ) j word2 ) * PITWOM31; 
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/* Rate of right ascension ( omega_dot ). 

*/ 

wordl = ( UNSIGNED TWOBYTE ) hex_ephemeris[22][0); 
MACRO_SET_UPPER(word1 1 ; /* Set upper bits if negative. */ 

MACRO_MAKE_EPH_WORD(word2, 22, 1, 22, 2 ); 

eph->omega_dot = ( ( ( FOURBYTE ) wordl << 16 ) | word2 ) * PITWOM43; 

I* Rate of inclination angle ( i_dot ). 

*1 

MACRO_MAKE_EPH_WORD( wordl, 23, 1, 23, 2 ) ; 
wordl = wordl > > 2 ; 

if ( wordl & 0x2000 ) /* Set upper bits if negative. */ 

wordl | = OxCOOO; 

eph->i_dot = ( TWOBYTE ) wordl * PITWOM43; 

I* Amplitude of the cosine harmonic correction term to the 
argument of latitude (cue). 

*/ 

M ACRO_M AKE_EPH_WORD(word 1 , 11,0, 11, 1 ); 
eph->cuc = ( TWOBYTE ) wordl * TWOM29; 

I* Amplitude of the sine harmonic correction term to the 
argument of latitude (cus). 

*/ 

M ACRO_M AKE_EPH_WORD(word 1 , 13, 0, 13, 1 ); 
eph->cus = ( TWOBYTE ) wordl * TWOM29; 

/* Amplitude of the cosine harmonic correction term to the 
orbit radius (crc). 

*/ 

MACRO_MAKE_EPH_WORD(word1 , 20, 0, 20, 1 ); 
eph->crc = ( TWOBYTE ) wordl * TWOM5; 

I* Amplitude of the sine harmonic correction term to the 
orbit radius (crs). 

*/ 

M ACRO_M AKE_EPH_WORD(word 1 , 8, 1, 8, 2 ); 
eph->crs = ( TWOBYTE ) wordl * TWOM5; 

I* Amplitude of the cosine harmonic correction term to the 
angle of inclination (cic). 

*/ 

MACRO_MAKE_EPH_WORD(word1 , 16, 0, 16, 1 ); 
eph->cic = ( TWOBYTE ) wordl * TWOM29; 

I* Amplitude of the sine harmonic correction term to the 
angle of inclination (cis). 

*/ 

MACRO_MAKE_EPH_WORD(word1 , 18, 0, 18, 1 ); 
eph->cis = ( TWOBYTE ) wordl * TWOM29; 
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r Estimated group delay differential ( tgd ). 

*1 

wordl = ( UNSIGNEDTWOBYTE ) hex_ephemeris[4][2]; 
MACRO_SET_UPPER(word1 ) ; /* Set upper bits if negative. */ 

eph->tgd = ( TWOBYTE ) wordl * TWOM31; 

I* Clock data reference time ( toe ). 

*/ 

MACRO_MAKE_EPH_WORD( wordl, 5, 1, 5, 2 ); 
eph->toc = wordl * 16.0; 

I* Polynomial coefficient ( af2 ). 

*/ 

wordl = ( UNSIGNED TWOBYTE ) hex_ephemeris[6][0]; 
MACRO_SET_UPPER( wordl ) ; /* Set upper bits if negative. */ 

eph->af2 = ( TWOBYTE ) wordl * TWOM55; 

I* Polynomial coefficient ( afl ). 

*/ 

MACRO_MAKE_EPH_WORD(word1 , 6, 1, 6, 2 ); 
eph->af1 = ( TWOBYTE ) wordl * TWOM43; 

I* Polynomial coefficient ( afO ). 

*/ 

wordl = ( UNSIGNED TWOBYTE ) hex_ephemeris[7][0]; 
MACRO_SET_UPPER( wordl ) ; /* Set upper bits if negative. */ 

MACR0_MAKE_EPH_W0RD(word2, 7, 1, 7, 2 ) ; 
word2 = word2 & OxFFFC ; /* clear 2 bits */ 

eph- > afO = ( ( I FOURBYTE ) wordl < < 16 ) | word2 ) * TWOM3125; 

/* Semi-major axis ( a ). 

*/ 

eph->a = eph->sqrt_a * eph->sqrt_a; 

I* Corrected mean motion ( n ). 

*/ 

eph->n = SQRMU / ( eph->a * eph->sqrt_a ) + eph->delta_n; 

I* Square root of eccentricity"2 subtracted from 1 ( esqrt ). 

*/ 

eph->e_sqrt = sqrt (1.0- eph->e * eph->e ); 

/* Omega dot - WE ( wkl ). 

*/ 

eph- > wkl = eph->omega_dot - WE; 

I* OmegaO - WE * toe ( wkO ). 

*/ 

eph->wkO = eph->omega_0 - WE * eph->toe; 

/* Sine w ( sin w ). 
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*/ 

eph->sin_w = sin ( eph->w ); 

I* Cosine w ( cos_w ). 

*/ 

eph->cos_w = cos ( eph->w ); 

I* Square root of semi-major axis * eccentricity * FCONST. 

*1 

eph->sqrta_e_f const = eph->sqrt_a * eph->e * FCONST; 

I* Issue of data, ephemeris. 

*/ 

eph->iode = hex_ephemeris[23)(0]; 

return *eph; 

}; 

/******♦*♦**♦* + + + *♦* + *** + ****♦*♦**♦♦♦♦ + ****♦**♦*♦******♦*******■•• 



Convert_Pos_Chan_Status() 

This procedure parses and decodes a position/status/data 
record received from the Motorola PVT-6. The flow of this 
procedure is as follows : 

- pass in the BINARY buffer containing 66 bytes 

- read the BINARY data from the buffer and convert directly 
into the final data structure. 

Input is a buffer of type ONEBYTE with a single binary 
position/status/data message. 

Output is the structure with the scaled and formatted 
position/status/data message. 

*1 

struct TPOSCHANSTATUS Convert_Pos_Chan_Status(ONEBYTE post]) 

{ 

struct TPOSCHANSTATUS pos_chan; 

UNSIGNED ONEBYTE i, j ; 

UNSIGNED_ONEBYTE tempchar ; 

UNSIGNED_FOURBYTE tempu4byte ; 

FOURBYTE temps4byte ; 

UNSlGNED_FOURBYTE nsecs ; 
double degrees, minutes, seconds ; 

/* Read and scale the of the data. 

*/ 

pos_chan. month = pos[2] ; 
pos_chan.day = pos[3] ; 

tempchar = pos(4] ; 
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pos_chan.year = ( tempchar < < 8 ) + pos[5] ; 

pos_chan. hours = pos[6] ; 
poschan. minutes = pos[7] ; 

tempchar = pos[8] ; /* integer seconds */ 
tempu4byte = pos[9] ; 

tempu4byte = ( tempu4byte << 8 ) + postIO] ; 

tempu4byte = ( tempu4byte < < 8 ) + pos[1 1 ] ; 

tempu4byte = ( tempu4byte < < 8 ) + pos[12] ; 

poschan. seconds = (double) tempchar + \ 

( ( (double) tempu4byte ) / 1.0E + 9 ); 

temps4byte = pos[13] ; 

temps4byte = ( temps4byte << 8 ) + post 14] ; 
temps4byte = ( temps4byte << 8 ) + pos[15] ; 

temps4byte = ( temps4byte < < 8 ) + pos[16] ; 

degrees = (double) temps4byte * MSECS_TO_DEGREES ; 

pos_chan. latitude. degrees = (TWOBYTE) degrees ; 

if ( degrees < 0 ) 

degrees = tabs ( degrees ) ; 

minutes = ( degrees - (TWOBYTE) degrees ) * 60.0 ; 
pos_chan. latitude. minutes = (TWOBYTE) ( minutes ) ; 
pos__chan. latitude. seconds = ( minutes - (TWOBYTE) minutes ) * 60.0 ; 

temps4byte = post 17] ; 

temps4byte = ( temps4byte << 8 ) + post 18] ; 

temps4byte = ( temps4byte < < 8 ) + post 19] ; 

temps4byte = ( temps4byte < < 8 ) + pos[20] ; 

degrees = (double) temps4byte * MSECS_TO_DEGREES ; 

pos_chan. longitude. degrees = (TWOBYTE) degrees ; 

if ( degrees < 0 ) 

degrees = fabs ( degrees ) ; 

minutes = ( degrees - (TWOBYTE) degrees ) * 60.0 ; 

pos_chan. longitude. minutes = (TWOBYTE) ( minutes ) ; 

pos_chan. longitude. seconds = ( minutes - (TWOBYTE) minutes ) * 60.0 

temps4byte = pos[21] ; 

temps4byte = ( temps4byte < < 8 ) + pos[22] ; 

temps4byte = ( temps4byte < < 8 ) + pos[23] ; 

temps4byte = ( temps4byte < < 8 ) + pos[24] ; 

pos_chan.datum_height = (double) temps4byte / 100.0 ; 

temps4byte = pos[25] ; 

temps4byte = ( temps4byte < < 8 ) + pos[26] ; 

temps4byte = ( temps4byte < < 8 ) + pos[27] ; 

temps4byte = ( temps4byte < < 8 ) + pos[28] ; 

pos_chan.msl_height = (double) temps4byte / 100.0 ; 
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tempchar = pos[29] ; 

pos_chan. velocity = (double) ( ( tempchar < < 8 ) + pos[30] ) / 100.0 ; 
tempchar = pos[31 ] ; 

pos_chan. heading = (double) ( ( tempchar < < 8 ) + pos[32] ) / 10.0 ; 
tempchar = pos[33] ; 

pos_chan.current_dop = (double) ( ( tempchar < < 8 ) + pos[34] ) / 10.0; 

pos_chan.dop_type = pos[35] ; 
pos_chan.visible_sats = pos[36] ; 
pos_chan.sats_tracked = pos[37] ; 

j = 0; 

for (i = 0; i < NUM_CHANNELS; i + + ) 

{ 

pos_chan.channel[i].svid = pos[38 + j]; 
pos_chan.channel[i].mode = pos[39 + j]; 
pos_chan.channel[i]. strength = pos[40 + j]; 
pos_chan.channel[i]. flags = pos[41 + j]; 
j+ =4; 

} 

pos_chan.rcvr_status = pos[621; 

return pos_chan; 

}; 

#endif /* _GPSC0NV_H */ 
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c 



GPSCRNCH • H 



I* GPSCRNCH. H 

Used to crunch the position data to determine the mean, 
maximum, and standard deviation of the position error 
obtained by DGPS. Geodetic coordinates are converted to 
ECEF coordinates and the distance formula is used to compute 
the distance error. 

Input is an array of structures that contain the scaled 
and formatted position/chan/status message data. 

The user is asked to supply the output file name 
without the extension. 

Output is a file with the distance error for each data 
point, and the mean, max and standard deviation of 
all the data points. 



tfifndef _GPSCRNCH_H 
^define J3PSCRNCH_H 

^include <math.h> 

^include <stdio.h> 

^include < string, h > 

#ifndef _GPSTRUCT_H 

^include n c:\borlandc\twite\gpstruct.h" 

#endif 

FILE *fp; 

void Ave_Max_Stdev(struct T_POS_CHAN_STATUS *mast, \ 

struct T_POS_CHAN_STATUS *slav, int master) 

{ 

int i, count = 0, start; 

char name[40], outpath[30] = ”c:\\borlandc\\twite\\"; 

double mlat, mlon, malt, slat, slon, salt, mN, sN, mean, distsum = 0; 

double mx, my, mz, sx, sy, sz, dist[200], max = -1.0, stdev, distsum2=0; 

printfC'Enter the output filename without the extension (???.txt): "); 

scanf("%s",name); 

strcatfoutpath, name); 

strcatfoutpath, ".out"); 

fp = fopen(outpath, M w M ); 

fprintf(fp, n \n%s\n",outpath); 

printfC'Enter the starting data point number: "); 
scanf("%d",&start); 



109 



for (i = start; i< master; i+ +) 

{ 

r Test to ensure the seconds match */ 

if ( (ONEBYTE) mast[i]. seconds = = (ONEBYTE) slav[i]. seconds ) 

{ 

r Convert to mx, my, mz, sx, sy, sz. 

*1 

mlat = (mast[i]. latitude. seconds * (1/60) + \ 

mast[i]. latitude. minutes) * (1/60) + mast[i]. latitude. degrees; 
mlon = (mastli]. longitude. seconds * (1/60) + \ 

mast[i]. longitude. minutes) * (1/60) + mastli]. longitude. degrees; 
malt = mast[i].datum_height; 

slat = (sla v[i] . latitude. seconds * (1/60) + \ 

slavli]. latitude. minutes) * (1/60) +slav[i]. latitude. degrees; 
slon = (slav[i] .longitude. seconds * (1/60) + \ 

slav[i]. longitude. minutes) * (1/60) + sla v[i] . longitude. degrees; 
salt = slav[i].datum_height; 

mN = SEMIMAJOR_AXIS / sqrt( 1 - E_SQUARED * sin(mlat) ); 
sN = SEMIMAJOR AXIS / sqrt( 1 - E_SQUARED * sin(slat) ); 

mx = ( mN + malt ) * cos(mlat) * cos(mlon); 
my = ( mN + malt ) * cos(mlat) * sin(mlon); 
mz = ( mN * ( 1 - E_SQUARED ) + malt) * sin(mlat); 

sx = ( sN + salt ) * cos(slat) * cos(slon); 

sy = ( sN + salt ) * cos(slat) * sin(slon); 

sz = ( sN * ( 1 - E_SQUARED ) + salt) * sin(slat); 

dist[i] = sqrt( (mx-sx)*(mx-sx) + (my-sy)*(my-sy) + (mz-sz)*(mz-sz) ); 

if ( distli] > 1000 ) 

{ 

fprintf(fp,"\nMtime = %2d:%10.7f <-> %2d:%10.7f = Stime \ 
Dist #%3d = XXXXX M ,mast[i]. minutes, mast[i]. seconds, \ 
slavli]. minutes, slav[i]. seconds, i); 
continue; 

}; 



I* Compute average, max and standard deviation. 

*/ 

count -i- + ; 

if ( distli] > max ) max = dist[i]; 
distsum + = distli]; 
distsum2 + = distli] * distli]; 

r Print out the results to a file. 

*/ 

fprintf(fp,"\nMtime = %2d:%10.7f <-> %2d:%10.7f = Stime \ 
Dist #%3d = %5.2f",mast[i]. minutes, mastli]. seconds, \ 
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slav[i]. minutes, slav[i]. seconds, i, dist[i]); 



else 

{ 

fprintf(fp,"\nMtime = %2d:%l0.7f <-> %2d:%10.7f = Stime \ 
Dist #%3d = XXXXX",mast[i]. minutes, mastli]. seconds, \ 
slavti]. minutes, slav[i] . seconds, i); 

}; 

}; 



/* Compute the mean and standard deviation. 

*/ 

mean = distsum / count; 

stdev = sqrt( (distsum2/count) - (mean * mean)); 

fprintfffp, "\n\nMean = %5.2f\nMax = %5.2f\nStdev = %5.2f\n",\ 
mean, max, stdev); 

fclose(fp); 

}; /* end Ave_Max_Stdev */ 

#endif /* GPSCRNCH H */ 
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D 



GPSDEFIN.H 



I* GPSDEFIN.H 



Store most of the definitions here. 

*/ 

tfifndef _GPSDEFIN_H 
^define _GPSDEFIN_H 

/********* type definitions *********/ 



typedef char ONEBYTE; /* data type allowing for the use of 8 bits. *1 
typedef short TWOBYTE; /* data type allowing for the use of 16 bits. */ 
typedef long FOURBYTE; /* data type allowing for the use of 32 bits. */ 

typedef unsigned char UNSIGNED_ONEBYTE; /* unsigned 8 bits */ 

typedef unsigned short UNSIGNED_TWOBYTE; /* unsigned 16 bits */ 

typedef unsigned long UNSIGNEDFOURBYTE; /* unsigned 32 bits */ 

/* ******** constant definitions *****/ 

^define PI 3.1415926535898 /* GPS value of Pi */ 

^define WE 7.2921 1 51 467E-5 /* WGS84 earth's rot rate (rads/sec) */ 

^define SQRMU 1 .9964981 84321 74E7 /* Square root of MU. */ 

^define FCONST -4.442807633E-1 0 /* FCONST is used in relativistic 

clock correction. Defined on 
page 73 of ICD-GPS-200. */ 

^define INVERSE_FLATTENING 298.257223563 /* WGS-84 unitless */ 

^define FLATTENING 1 /INVERSEFLATTENING 

^define SEMIMAJOR_AXIS 63781 37.0 /* WGS-84 meters */ 

^define E SQUARED 2* FLATTENING-FLATTENING ‘FLATTENING 



^define TWOM5 


0.03125000000000 


^define TWO Ml 1 


(TWOM5 / 


64.0) 


^define TWOM19 


(TWOM1 1 / 256.0) 


#define TW0M20 


(TWOM19 / 


2.0) 


^define TWOM21 


(TWOM20 / 


2.0) 


^define TWOM23 


(TWOM21 / 


4.0) 


^define TWOM24 


(TWOM23 / 


2.0) 


^define TWOM27 


(TWOM24 / 


8.0) 


^define TWOM29 


(TWOM23 / 


64.0) 


^define TWOM30 


(TWOM29 / 


2.0) 


^define TWOM31 


(TWOM30 / 


2.0) 


^define TWOM33 


(TWOM31 / 


4.0) 


^define TWOM38 


(TWOM33 / 


32.0) 


^define TWOM43 


(TWOM38 / 


32.0) 
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^define TWOM50 (TWOM43 / 1 28.0) 

^define TWOM55 (TWOM50 / 32.0) 

^define TWOM3125 (0.25 * TWOM31) 



^define TWO P1 1 
^define TW0P12 
^define TW0P14 
^define TW0P1 6 



2048.0 /* 2*11 */ 

(TW0P1 1 * 2.0) 

(TW0P12 * 4.0) 

(TW0P14 * 4.0) 



//define TW0P1 2JNTEGER 4096 /* 2" 12 */ 

//define TW0P1 6JNTEGER 65536 /*2*16 */ 

//define PITW0M31 (PI * TW0M31) 

//define PITW0M43 (PI * TWOM43) 



/********* Twite definitions *♦****♦*/ 



//define 

//define 

//define 

//define 

//define 

//define 

//define 

//define 



FALSE 
TRUE 1 
RX 0 

TX 1 

MGPS_PORT 
SGPSPORT 
NUM_CHANNELS 6 
MSECS TO DEGREES 



/* receive */ 
I* send */ 

3 

4 



I* master gps receiver port *1 



I* slave gps receiver port */ 
((1.0 / 1000.0 ) / 3600.0 ) 



//define EPH_NUM_WORDS 24 /* at 3 bytes/word */ 

//endif /* GPSDEFIN H */ 
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E 



GPSDELAY.H 



/* GPSDELAY.H 

Used to introduce an artificial delay in transmitting 
psuedorange corrections from the master to the slave. 

PR_delay() accepts a pointer of type PANDL to the buffer 
containing the psuedorange correction message 
(@@Ce...) and an integer number of seconds it 
is to be delayed. The output is to the PCL-744 
serial I/O card. 

*/ 

#ifndef _GPSDELAY_H 
^define _GPSDELAY_H 

include <stdio.h> 

#include <string.h> 

^include < stdlib. h > 

^include <time.h> 

^include <dos.h> 

^include "c:\borlandc\twite\head-c.h" 

#ifndef _GPSTRUCT_H 

^include "c:\borlandc\twite\gpstruct.h" 

#endif 

^define MAXDELAY 300 

static PANDL stor[MAXDELAY]; /* Ring buffer */ 

static int num_waiting = 0; /* Number of messages in queue */ 

static int next_stor = 0; /* Next storage location to be used */ 
static int next_xmit = 0; I* Next message to be transmitted */ 

void PR_delay( PANDL *Ce_msg, int delay ) 

{ 

PANDL *ptr; 

I* Copy the contents of the buffer with the Ce msg to a 
new buffer for holding during the delay. 

*/ 

ptr = mallocl sizeof(PANDL) ); 
ptr->len = Ce_msg->len; 

ptr->ptr = callocf ptr->len, sizeof(ONEBYTE) ); 
memcpyf ptr->ptr, Ce_msg->ptr, Ce_msg->len ); 

/* Store the structure of type PANDL in the ring buffer ”stor[]". 

*/ 

if (num waiting < MAXDELAY ) 
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{ 

stor[next_stor] = *ptr; 
numvvaiting + + ; 
next_stor + + ; 

if (nextstor = = MAXDELAY) 
next_stor-= MAXDELAY; 

} 

else 

{ 

printfC'There are %d \"@@Ce...\” messages in the ring buffer. \n",\ 
MAXDELAY); 

printf("One correction message was lost.”); 

free(ptr->ptr); 

free(ptr); 

return; 

} 

r Transmit the oldest message if there are more messages 
than the delay length. The delay is premised upon the 
fact that correction messages come approximately one second 
apart. 

*/ 

while ( (num_waiting - delay) > 0 ) 

{ 

sio_putb( SGPS_PORT, stor[next_xmit].ptr, stor(next_xmit].len ); 

free(stor[next_xmit].ptr); 

next_xmit + + ; 

if (next_xmit = = MAXDELAY) 
next_xmit - = MAXDELAY; 
num_waiting- ; 

}; 



free(ptr); 

return; 

}; 

#endif /* _GPSDELAY_H */ 
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F 



GPSFUN.H 



/* GPSFUN.H 

A collection of functions to manipulate GPS messages. 
Hex() 

Check_sum() 

Bin^oasciiO 

Msg_verification() 

ParsejnbuffO 

Mastergpsf) 

Slave_gps() 

*/ 

#ifndef _GPSFUN_H 
#define _GPSFUN_H 

^include < stdlib.h > 

^include <string.h > 

^include <stdio.h> 

^include "c:\borlandc\twite\head-c.h" 

#ifndef _GPSTRUCT_H 

^include "c:\borlandc\twite\gpstruct.h" 

#endif 

FILE *m_all_file; /* master file for all data */ 

FILE *m_pos_file; /* master file for pos data */ 

I* Global matrix for converting ephemeris 
*1 

UNSIGNED_ONEBYTE hex_ephemeris[24][3] ; 

/* [word][byte] : 24 words (3 bytes per word) corresponding to 
words 3-10 of subframes 1-3 of the satellite nav message. 

*/ 



yf*** + *** + ******** + + *** + * + + ***** + + + + + * + + * + + + + + + + + + + * + + + **** + * + + * + 



Hex 

This procedure converts an ASCII value to its hex 
equivalent. 

Input is a single ONEBYTE. 

Output is a single ONEBYTE. 
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ONEBYTE HexlONEBYTE ASCIIvalue) 

{ 

/* 0 is 0x30(48) -> 9 is 0x39(57) 

A is 0x41(65) -> F is 0x46(70) 
a is 0x61(97) -> f is 0x66(102) 

*/ 

if ((ASCIIvalue) < 0x40) 

return((ASCIIvalue) - 0x30); 

else 

if ((ASCIIvalue) < 0x60) 

return((ASCIIvalue) - 0x37); 

else 

return((ASCIIvalue) - 0x57); 

}; 



y***************‘ 



Check_sum() 

Calculates the checksum of a message in a buffer of type 
ONEBYTE. Since the message is assumed to be in the Motorola 
Binary Proprietary format, all but the last three bytes 
are used in computing the checksum. 

Input is the pointer of type PANDL. 

Output is the checksum of type ONEBYTE. 



ONEBYTE Check_sum( PANDL *msg ) 

{ 

ONEBYTE chksum = 0; 
int i; 

I* Compute the Exclusive-Or of all but the last three bytes 
*1 

for ( i = 0; i < (msg->len - 3); i+ + ) 
chk_sum " = msg->ptr[i]; 
return chk_sum; 

}; 



/ 



**+**#*******♦♦***+****++*++++*++***+******+******+**+**♦♦+♦*+* 



Bin_to_ascii() 

Converts a message from binary to its hexidecimal 
ASCII equivalent. All bytes of the input message are 
converted to ASCII except the first n values as specified 
in the integer input argument. (This was incorporated to 
permit the use of this function with the Motorola GPS 
proprietary message format where the first four bytes are 
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already in ASCII.) 



Input is a pointer of type PANDL and the interger value 
of the number of bytes to skip before conversion. 
Output is a pointer of type PANDL to the new buffer 
with the ASCII converted message. 



PANDL *Bin_to_ascii( PANDL *bin_msg, int skip ) 

{ 

int i, j, index, tmplen = 0 ; 

ONEBYTE nibble[2], *ptr_tool; 

PANDL *ascii_msg; 

tmplen = bin_msg->len * 2 - skip; 
ptr_tool = bin_msg->ptr; 

r Create new memory for the converted message. 

*/ 

ascii_msg = mallocf sizeoff PANDL )); 
ascii_msg->ptr = calloc( tmplen, sizeof(ONEBYTE) ); 
ascii_msg-> len = tmplen; 

r Copy any initial values to the new buffer that don't need conversion. 
*/ 

memcpy(ascii_msg->ptr,ptr_tool,skip); 

for ( i = skip; i < tmplen; i + + ) 

{ 

/* Split the byte into its upper and lower nibbles. 

V 

nibbletO] = ( ptr_tool[i] > > 4 ) & OxOF; /* upper */ 

nibbled 1 = ( ptr_tool[i] & OxOF ); /* lower */ 

for ( j = 0; j < 2; j + + ) 

{ 

I* Convert the nibbles to ASCII. 

*/ 

if ( nibblelj] < 10 ) 

nibblelj] + = 0x30; /* nibble + 0x30 = 0x30:0x39 */ 

else 

nibblelj] += 0x37; /* nibble + 0x37 = 0x41 :0x46 */ 

}; 

r Compute the index and insert into the ascii_msg buffer 
*/ 

index = 2 * i - skip; 

ascii_msg->ptr[ index ] = nibblelO]; 

ascii_msg->ptr[ index + 1] = nibble! 1 ] ; 

}; 



return ascii_msg; 
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Msg_verification() 

Specifically used to verify the integrity of the 

Motorola GPS proprietary message format. 

Input are a pointer of type PANDL and an integer designating 
whether the message was received from the GPS receiver (Rx) 
or is to be transmitted to the GPS receiver (Tx). 

Output is the integer flag indicating whether or not the 
message passed verification. True means it passed. 



int Msg_verification( PANDL *msg, int msg_direction ) 

{ 

int i, len = -1 , arraycol; 

int status = TRUE; 



ONEBYTE 

ONEBYTE 

ONEBYTE 

ONEBYTE 

ONEBYTE 

ONEBYTE 

ONEBYTE 



byteO = msg->ptr[0]; /* Msg prefix @ byte */ 

bytel = msg->ptr[1]; /* Msg prefix @ byte */ 

byte2 = msg->ptr[2]; /* Msg id first byte */ 

byte3 = msg->ptr[3]; /* Msg id second byte */ 

byteCS = msg->ptr[ msg->len-3 ]; /* Check Sum byte */ 
byteCR = msg->ptr[ msg->len-2 ]; /* Carr Retn byte *7 
byteLF = msg->ptr[ msg->len-1 ]; /* Line feed byte */ 



I* Determine if message verification is for Tx or Rx because the 
number of byte in the message is different for a Tx message 
than for an Rx message. 

*/ 

if ( msgdirection ) 

array_col = 0; /* Tx message: to GPS receiver */ 

else 

array col = 3; /* Rx message: fm GPS receiver */ 

I* Message ID existence verification and expected length. 

*/ 

if ( byte2 = = 'A' ) 

for ( i = 0; i < 30; i + + ) 

{ 

if ( byte3 = = A_array[i][array_col] ) 

{ 

if ( msg direction ) 

len = A_array[i][1 ]; l* Msg size to PVT-6 */ 
else 

len = A_array[i][2]; /* Msg size fm PVT-6 */ 
break; 

}; 
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}; 

if ( byte2 = = 'B' ) 

for (i = 0; i < 1 1 ; i + + ) 

{ 

if ( byte3 = = B_array[i][array_col] ) 

{ 

if ( msgdirection ) 

len = B_array[i][1 ]; /* Msg size to PVT-6 *1 

else 

len = B_array[i][2]; /* Msg size fm PVT-6 */ 
break; 

}; 

}; 

if ( byte2 = = 'C' ) 

for ( i = 0; i < 1 1 ; i + + ) 

{ 

if ( byte3 = = C_array[i][array_col] ) 

{ 

if ( msg direction ) 

len = C_array[i][1 ]; I* Msg size to PVT-6 */ 
else 

len = C_array[i][2]; /* Msg size fm PVT-6 *1 
break; 

}; 

}; 

if ( byte2 ! = 'A' && byte2 ! = 'B' && byte2 ! = X' ) 

{ 

status = FALSE; 
return status; 

}; 



/* Message length verification. (”Cj” message is a special case. 
See GPSINIT.H.) 

*/ 

if (len != msg->len && byte2 != X' && byte3 != 'j' ) 

{ 

status = FALSE; 
return status; 

}; 

/* Message check sum verification. 

*1 

if ( byteCS I = Check sum! msg ) ) 

{ 

status = FALSE; 
return status; 

}; 

I* Message begins with @@ prefix. 

*/ 

if ( byteO !='<§>' j| bytel ! = '@' ) 
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{ 

status = FALSE; 
return status; 

}; 

r Message ends with CR-LF. 

*/ 

if ( byteCR ! = OxOD ] | byteLF ! = OxOA ) 

{ 

status = FALSE; 
return status; 

}: 

return status; 

}: 






Parse_inbuff() 

This function accepts a pointer of type PANDL to the 
contents extracted from the serial port Rx buffer. As 
this function repeatedly parses each succeeding message 
from this buffer, the buffer's length decreases until 
the only remaining information is a partial message. 

The pointer is a static pointer and will not lose its 
information upon completion of this function so it will 
be rejoined with the rest of itself when the next grab 
from the Rx buffer occurs. 

Inputs are a pointer of type PANDL and a pointer to a flag 
indicating when the only remaining contents constitute 
a partial message. 

Output is a pointer to a new buffer of type PANDL containing 
the next complete single message in the buffer. Also, 
the inbuff pointer passed into Parsejnbuff has been 
altered. 



PANDL *Parse_inbuff( PANDL *start_of_next_message, int *partial_msg_flag ) 

{ 

int i = 0, j = 0; 

int msg_start_found_before = FALSE; 

PANDL *parsed_msg; 

ONEBYTE *ptr_tool; 

parsed msg = mallocf sizeoft PANDL ) ); 

ptr_tool = start_of_next_message->ptr; 

r CASE #1: A complete message has been found. The pointer 
tool points to the beginning of a msg (indicated by 
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and "@@" has been found once before. Move the place keeper 
pointer ”start_of_next_msg->ptr H and adjust the length 
”start_of_next_msg->len". Copy the "found" message to a 
new buffer and return its location in a pointer of type 
PANDL. 

*/ 

for ( i = 0; i < start_of_next_message->len; i + + ) 

{ 

if ( ptr_tool[i] = = && ptr_tool[i + 1] = = ) 

{ 

if ( msgstartfoundbefore ) 

{ 

parsed_msg->len = i - j; 

parsed_msg->ptr = calloc(parsed_msg-> len, sizeof(ONEBYTE)); 
memcpy(parsedjnsg->ptr, ptrtool, parsed jnsg-> len); 

start_of_next_message->ptr + = i; 
start_of_next_message->len - = i; 

*partialjnsg_flag = FALSE; 
return parsedjnsg; 

}; 

msg_start_found_before = TRUE; 
j = i; 

}; 

}; 

r The start of a second message "@@" was not found. The contents 
are either exactly a complete message or they constitute a 
partial message. Thus, the Parsejnbuff function is 
given another chance otherwise the "partial-msg_flag" is 
returned and the partial message contents will be prepended 
to the next buffer brought into Parsejnbuff. 

*/ 

r CASE #2\ There is exactly one message in the buffer. 

*/ 

if ( ptr tooKO] = = '@' && \ 
ptr tool! 1 1 == && \ 

ptr_tool[ startjDfjiextjriessage->len - 2 ] = = OxOD && \ 
ptr_tool[ start_of_next_message-> len - 1 ] = = OxOA ) 

{ 

parsed_msg->len = i - j; 

parsed_msg-> ptr = calloc(parsed_msg-> len, sizeof(ONEBYTE)); 
memcpy(parsed_msg->ptr, ptr_tool, parsedjnsg->len); 

start_of_next_message-> ptr = NULL; 
startjDfjiext_message->len = 0; 

*partial_msg_flag = FALSE; 

} 
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/* CASE #3: There is only a partial message left in the buffer. 
*/ 
else 

{ 

*partial_msg_flag = TRUE; 

}; 

return parsed_msg; 

}; 



/ 



+*+****++**++**+*+***+♦*+*+*++*++**********++*++++*+*++***++*** 



MasterjgpsO 

For use with a Motorola GPS receiver configured in the 
master station mode. It builds two binary files on disk: 

One records all output messages and the other records just 
position messages. The function returns a pointer of type 
PANDL that points to a buffer with psuedorange correction 
messages ready for transmission to the slave receiver. 

Input is nothing. 

Output is a pointer of type PANDL. 



PANDL *Master_gps(void) 

{ 

static PANDL mxsbuff; 

static ONEBYTE mxs[3001; /* buffer for partial msg *1 

static ONEBYTE mstart = TRUE ; 

ONEBYTE * tail; 

PANDL *inbuff, *tmpbuff, *parse_tool, *ben; 

int *partial_msg_flag; 

partial_msg_flag = mallocl sizeofl ONEBYTE ) ); 
*partial_msg_flag = TRUE; 

inbuff = mallocl sizeofl PANDL ) ); 
tmpbuff = mallocl sizeofl PANDL ) ); 
parse_tool = mallocl sizeofl PANDL ) ); 
ben = mallocl sizeofl PANDL ) ); 

I* Determine amount of info in Rx buffer and import to tmpbuff. 
*1 

tmpbuff- >len = (ONEBYTE) siojqueuel MGPS PORT ); 
tmpbuff- >ptr = callocl tmpbuff- >len, sizeof(ONEBYTE) ); 
sio_read( MGPS_PORT, tmpbuff- >ptr, tmpbuff- >len ); 

/* Initialize the excess buffer that is to contain the partial msg. 
*/ 

if ( mstart ) 

{ /* This "if” statement is invoked *1 
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mxsbuff.len = 0; /* only the first time the program */ 

mstart = FALSE; /* is run. Note the scope of mstart.*/ 
mxsbuff.ptr = mxs; 

}; 

f* Append tmpbuff messages to the partial message that is in mxsbuff. 
*/ 

inbuff- >len = tmpbuff- >len + mxsbuff.len; 
inbuff->ptr = callocl inbuff->len, sizeof(ONEBYTE) ); 
memcpyl inbuff- >ptr, mxsbuff.ptr, mxsbuff.len ); 
tail = inbuff->ptr + mxsbuff.len; 
memcpyl tail, tmpbuff->ptr, tmpbuff->len ); 

r Free tmpbuff memory and reset partial message buffer length to zero. 
*/ 

free! tmpbuff- >ptr ); 
free! tmpbuff ); 
mxsbuff.ptr = mxs; 
mxsbuff.len = 0; 

I* Initialize "ben" for use in collecting psuedorange correction messages. 
*/ 

ben = inbuff; 
ben->len = 0; 

while ( TRUE ) 

{ 

I* Get a full message from inbuff. 

*1 

parse tool = Parse inbuffl inbuff, partial msg flag ); 
if ( *partial_msg_flag ) 

{ 

free! parse_tool->ptr (; 
break; 

}; 



I* Verify the received message validity and indicate when 
the message fails verification. 

*1 

if ( Msg_verification( parse_tool, RX ) ! = TRUE ) 

{ 

printfC'Message from master receiver fails verification!\n"); 

freel parse_tool->ptr ); 

continue; 

}; 



I* Store only psuedorange corrections in old inbuff buffer for 
return to the calling routine. 

*/ 

if (parse_tool-> ptr[2] = = 'C' && parse_tool->ptr(3] = = 'e' ) 

memcpy(ben->ptr + ben->len, parse_tool-> ptr, parse tool- > len); 
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I* Write all messages to this file. 

*/ 

m_all_file = fopen( "mall.txt", "ab" ); 

fwritefparsetool- > ptr,sizeof(ONEBYTE),parse_tool- > len,m_all _file); 
fclosef m_all_file ); 

/* Write only position messages to this file. 

*1 

if (parse_tool- > ptr[2] = = 'B' && parse_tool->ptr[3] = = 'a' ) 

{ 

m_pos_fil e = fopenf "mpos.txt", "ab" ); 

fwrite(parse_tool- > ptr,sizeof(ONEBYTE),parse_tool- >len,m_pos_file); 
fclose( m_pos_file ); 

}; 



free( parse_tool->ptr ); 

}; /* end while */ 

I* Store any partial message in the excess buffer. 

*/ 

if ( inbuff- >len > 0 ) 

{ 

memcpyfmxsbuff.ptr, inbuff->ptr, inbuff- >len); 
mxsbuff.len = inbuff->len; 

}; 



free( inbuff ); 
free! parse_tool ); 
free( partial msg flag ); 

return ben; 

}; 






Slave_gps() 

For use with a Motorola GPS receiver in the slave 
station mode. It converts five data messages into 
scaled and formatted data: 

Position/Status/Data message 
xDOP Table Status message 
Ephemeris Data Output message 
Satellite Range Data Output message 
Position/Status/Data Extension message 
In addition, these and all other output messages are 
returned to the calling routine for transmission to 
the master receiver for post-processing. 

Input is a pointer of type PANDL with any input messages. 
Output is a pointer of type PANDL with all output messages. 



125 



/ 



PANDL *Slave_gps( PANDL *input) 

{ 

static PANDL sxsbuff; 

static ONEBYTE sxs[300]; I* buffer for partial msg */ 
static ONEBYTE sstart = TRUE ; 

ONEBYTE *tail; 

PANDL *inbuff, *tmpbuff, *parse_tool, *ben; 

int *partial_msg_flag, i; 

r Deliver any input messages to the GPS receiver. 

*/ 

sio_putb( SGPS_PORT / input->ptr, input->len ); 
free( input->ptr ); 
free( input ); 

r Create flag pointer and working pointers of type PANDL. 

*/ 

partial_msg_flag = malloc( sizeof(ONEBYTE) ); 

*partial_msg_flag = TRUE; /* Used in Parsejnbuff */ 

inbuff = malloc( sizeoft PANDL ) ); 
tmpbuff = malloct sizeof( PANDL ) ); 
parse_tool = malloct sizeoft PANDL ) ); 
ben = malloct sizeoft PANDL ) ); 

/* Determine amount of info in Rx buffer and import to tmpbuff. 

*/ 

tmpbuff- >len = (ONEBYTE) sio_iqueue( SGPS PORT ); 
tmpbuff->ptr = calloct tmpbuff-> len, sizeof(ONEBYTE) ); 
sio_read( SGPS_PORT, tmpbuff->ptr, tmpbuff->len ); 

I* Initialize the excess buffer that is to contain the partial msg. 

*1 

if ( sstart ) 

{ r This "if" statement is invoked */ 

sxsbuff. len = 0; I* only the first time the program */ 
sstart = FALSE; /* is run. Note the scope of sstart. */ 
sxsbuff. ptr = sxs; 

}; 



r Append this info in tmpbuff to excess msg from last read. 

* / 

inbuff->len = tmpbuff-> len + sxsbuff. len; 
inbuff->ptr = calloct inbuff->len, sizeof(ONEBYTE) ); 
memcpyt inbuff- > ptr, sxsbuff. ptr, sxsbuff. len ); 
tail = inbuff->ptr + sxsbuff. len; 
memcpyt tail, tmpbuff-> ptr, tmpbuff->len ); 

r Free tmpbuff memory and reset partial message buffer length to zero. 



126 



*1 

f ree< tmpbuff->ptr ); 
free( tmpbuff ); 
sxsbuff.ptr = sxs; 
sxsbuff.len = 0; 

I* Initialize "ben" for use in returning all messages to the calling 
environment. 

V 

ben = inbuff; 

while ( TRUE ) 

{ 

I* Get a full message from inbuff. 

*/ 

parse_tool = Parse_inbuff( inbuff, partial_msg_flag ); 
if ( *partial_msg_flag ) 

{ 

free! parse_tool->ptr ); 
break; 

}; 

I* Verify the received message validity. 

*/ 

if ( Msgjverificationl parse_tool, RX ) ! = TRUE ) 

{ 

free! parse_tool- > ptr ); 
continue; 

}; 



/* Convert the binary msg to GPS data structure for use by cpu. 
*/ 

switch (parse_tool-> ptr(31) 

{ 

case 'a'; /* @@Ba Position/Status/Data */ 

/* Insert conversion here */ 
break; 

case 'c': /* @@Bc xDOP Table Status */ 

I* Insert conversion here */ 
break; 

case 'f': /* @@Bf Ephemeris Data */ 

for ( i = 5 ; i < EPH_NUM_WORDS + 5 ; i+ =3 ) 

{ 

hex_ephemeris[i][0] = parse_tool-> ptr[i + OJ ; 
hex_ephemeris[i][1 ] = parse_tool->ptr[i + 1] ; 
hex_ephemeris[i][2] = parse_tool->ptr[i + 2] ; 

} 

r Insert conversion here */ 
break; 

case 'g': /* @@Bg Satellite Range */ 

!* Insert conversion here * I 
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break; 

case 'k': f* @@Bk Position/Status/Data Ext */ 

r Insert conversion here */ 
break; 
default: 
break; 



free( parse_tool->ptr ); 

}; /* end while */ 

/* Store any partial message in the excess buffer. 

*/ 

if ( inbuff- >len > 0 ) 

{ 

memcpy(sxsbuff.ptr, inbuff- >ptr, inbuff- > len); 
sxsbuff.len = inbuff->len; 

}; 

free( inbuff ); 
free( parse_tool ); 
free( partial_msg_flag ); 

return ben; 

}; 

tfendif /* GPSFUN H */ 
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G 



GPSINIT • H 



I* GPSINIT. H 

This header file contains the Motorola Proprietary 
Binary Message Format characterises for GPS 
receivers. The first column contains the second 
letter of the message ID of a message prepared for 
input to the GPS receiver. The fourth column contains 
the second letter of the message ID of a message 
received from the GPS receiver. For all messages sent 
to the GPS receiver except one, there is a message reply. 

The second column contains the final length of a 
message prepared to be sent to the GPS receiver. The 
third column is the length of a message received from 
the GPS receiver with that particular message ID. 

Note: All message lengths can be expressed with 
one byte except one: The response message Cj is 294 bytes. 
A unique dummy size equal to 1 byte has been used to 
indicate the length. 

*1 

tfifndef _GPSINIT_H 
^define _GPSINITJH 

tfifndef _GPSDEFIN_H 

^include ,, c:\^?o^landc\twite\gpsdefin.h ,, 

tfendif 



ONEBYTE 



A_array[30][4] = {{ 


'a' 


{ 


'b', 


, 10, 


10, 


'b' }, 


{ 


'c', 


,11, 


11, 


'c' }, 


{ 


'd', 


, 11, 


11, 


'd' }, 


{ 


'e', 


,11, 


11, 


'e' }, 


{ 


'f\ 


12, 


15, 


'f' }- 


{ 


'g', 


, 8, 


8, 


'g' }. 


{ 


'h', 


, 8, 


8, 


'h' }, 


{ 


/ j » 


9, 


9, ' 


}, 


{ 


r! 


8, 


8, ' 


j' }, 


{ 


'k\ 


, 9, 


9, 


'k' }. 


{ 


T, 


9, 


9, ' 


1' }. 


{ 


'm' 


’, 12 


- 12 


, 'm' } 


{ 


'n' 


, 8, 


8, 


'n' }, 


{ 


'o' 


, 8, 


25, 


'P' }, 


{ 


'p' 


, 25, 


25, 


'P' }. 


{ 


V 


, 8, 


8, 


'q' }. 


{ 


'r', 


8, 


8, ' 


’r' }, 


{ 


's', 


, 20, 


20, 


's' }, 


{ 


't', 


8, 


8, 1 


’t' }, 
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ONEBYTE 



ONEBYTE 



#endif I* 



{ 


V, 


12, 


12, 


'u 


' }, 






{ 


V, 


8, 


8, 1 


V 


}. 






{ 


'w' 


, 8, 


8, 


'w 


' }. 






{ 


'x', 


9, 


9, ‘ 


x' 


}. 






{ 


V, 


11. 


11, 


'y 


' }. 






{ 


’z\ 


11. 


11, 


'z' 


}. 






{ 


'A' 


8, 


8, 


'A' 


}, 






{ 


'B' 


8, 


8, 


'B' 


). 






{ 


'C' 


9, 


9, 


'C' 


). 






{ 


'D' 


9, 


9, 


'D' 


}}; 






_arrayl11H4] = {{ 




'a' 


od 

CD 

od 


'a' }, 


{ 


'b' 


8, 


92, 


'b' 


}, 






{ 


'c', 


8, 


82, 


'c' 


}- 






{ 


'd' 


8, 


23, 


'd' 


}. 






{ 


'e\ 


8, 


0, ' 




}./* 


Be->Cb * 


/ 


{ 


'i'. 


8, 


80, ' 


f' 


}. 






{ 


'9', 


, 8, 


122, 


'g' 


}. 






{ 


'h', 


, 8, 


0, ■ 




}./* 


Bh->Ce * 


7 


{ 


T, 


7, 


8, 'j 


r }, 






{ 


T, 


80, 


0, 




}./* 


Bf-> Cc * 


/ 


{ 


'k', 


8, 


69, 


'k' 


}}; 






:_array[1 1][4] = {{ 


f * 


Y, 7, 


9, 'a' 


}. 


{ 


'b', 


, 33, 


9, 


'h' 


}, 






{ 


'd', 


, 27, 


171 


,'d 


' }. 






{ 


'e', 


52, 


7, 


'k' 


}. 






{ 


r, 


7, 


7, ' 


V 


}, 






{ 


'g', 


, 8, 


8, 1 


’g' 


}. 






{ 


7, 


8, 


0, 'i 


i' ] 


f. 






{ 


7, 


7, 


1, 'j 


i' ] 


f, /* 


Cj resp msg 294 byte */ 


{ 


$ t 


0, 


33, 


'b' 


}< / 


* Be- > Cb 


*1 


{ 




0, 


52, 


'e' 


}. /■ 


* Bh->Ce 


*/ 


{ 


# t 


0, 


80, 


'c' 


}};/ 


* Bf->Cc 


*/ 



GPSINIT_H */ 
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H. 



GPSPORTS . H 



r GPSPORTS. H 

lnitialize_master_gps() 

lnitialize_slave_gps() 

OpenportsO 

Close_ports() 

Both initialization functions need to be greatly improved 
to provide an easy to understand and use interface with 
the user. 

*/ 

tfifndef _GPSPORTS_H 
^define _GPSPORTS_H 

include <stdio.h> 

^include "c:\borlandc\twite\head-c.h" 

r ********************************************** 



lnitialize_master_gps() 

This function presents the current configuration 
settings for the master GPS receiver. The user is 
offered the opportunity to change any of the parameters. 
In addition, the user is requested to indicate which if 
any output files he will want. For those files, the default 
file name is presented and made available for modification 
or change. 

Input is the serial port number through which the master 
GPS communicates. 

Output is nothing. 

*/ 

void lnitialize_master_gps(int mgps) 

{ 

char si [20] = {'@y@YBYa',1 ,22,13,10}; 
char s2[20] = {'@Y@YB','h',1 ,43, 1 3,1 0}; 
sio_flush{mgps,2); 
sio_putb(mgps,s1 ,8); 
sio_putb(mgps,s2,8); 

printf("Master GPS receiver initialized. \n"); 
return; 

}; 



+ * + + + + ♦* + + * + + + * + + * + + + + ♦ + ***♦** 
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Initialize_slave_gps() 

This function presents the current configuration 
settings for the slave GPS receiver. The user is 
offered the opportunity to change any of the parameters. 
In addition, the user is requested to indicate which if 
any output files he will want. For those files, the default 
file name is presented and made available for modification 
or change. 

Input is the serial port number through which the slave 
GPS communicates. 

Output is nothing. 



*/ 



void lnitialize_slave_gps(int sgps) 

{ 

char si [20] = ,22,1 3,1 0}; 

sio_flush(sgps,2); 

sio_putb(sgps,s1 ,8); 

printfC’Slave GPS receiver initialized. \n"); 
return; 

}; 



I * 1 



I- ****** H 



Open_ports() 

This function initializes the two serial ports 
on the PCD-744 to be used by the master and slave GPS 
receivers. Before initialization, it either loads the 
drivers for the standard and PCD-744 serial ports or 
verifies that they are loaded. It then configures the 
Baud rating, word length, parity, and stop bits as well 
as the modem and flow control. All commands that begin 
with ,, sio_...” are PCD-744 commands and may be referenced 
in the PCLS-802 PC-ComLIB Serial Communication Programming 
Library. 

The PCD-744 RS-232 8 port intelligent 10 card is 
made by ADVANTECH in Sunnyvale, Ca. (408) 245-6678. 

Input is the slave GPS port number. 

Output is nothing. 



void Open_ports( int mgps, int sgps ) 

{ 

int BAUD =12; /* 9600 (p. 1 2 PC-ComLIB) */ 

int IOMODE = (0x03 1 0x00 1 0x00); /* 8-N-1 (p. 12) */ 

int MODMODE = 0x00; /* DTR and RTS off (p.26)*/ 
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/' HW and SW flow Ctrl off 



int HWMODE 
int n; 



= 0x00; 
(p.33)*/ 



n = sio_ioctl( mgps, BAUD, IOMODE ); 
if (n ! = 0) 

{ 

printf("Opening comport %d ioctl error. \n", mgps); 
abortO; 

}; 

n = sio_lctrl( mgps, MODMODE ); 
if (n ! = 0) 

{ 

printff Opening comport %d Ictrl error. \n", mgps); 
abortO; 

}; 

n = sio_flowctrl( mgps, HWMODE ); 
if (n != 0) 

{ 

printf("Opening comport %d flowctrl error. \n” # mgps); 
abortO; 

}; 

n = sio_open( mgps ); 
if (n != 0) 

{ 

printfC'Opening comport %d error. \n”, mgps); 
abortO; 

}; 



n = sio_ioctl( sgps, BAUD, IOMODE ); 
if (n ! = 0) 

{ 

printfC'Opening comport %d ioctl error.\n",sgps); 
abortO; 

}; 

n = sio_lctrl( sgps, MODMODE ); 
if (n != 0) 

{ 

printfC'Opening comport %d Ictrl error. \n", sgps); 
abortO; 

}; 

n = sio_flowctrl( sgps, HWMODE ); 
if (n != 0) 

{ 

printfC'Opening comport %d flowctrl error. \n", sgps); 
abortO; 

}; 

n = sio_open( sgps ); 
if (n ! = 0) 

{ 
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printf ( ” Opening comport %d error. \n", sgps); 
abort!); 

}; 

printfC'Opened ports %d (Master GPS) and %d (Slave GPS).\n", 

mgps, sgps); 



return; 

}; 



\ 



Close_ports() 

This function closes and resets all serial ports for 

the PCD-744 Serial 10 card. See the Open_ports function 

for more information on the commands used herein. 

Inputs are the master and slave GPS port numbers. 
Output is nothing. 

*/ 

void Close_ports( int mgps, int sgps ) 

{ 

printfCCIosing all ports\n"); 
sio_close( sgps ); 
sio_close( mgps ); 
sio_reset(); 
return; 

}; 

tfendif /* _GPSPORTS_H */ 
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GPSTRUCT • H 



r GPSTRUCT. H 
*/ 

tfifndef _GPSTRUCT_H 
^define _GPSTRUCT_H 

tfifndef _GPSINIT_H 

^include " c:\borlandc\twite\gpsinit. h" 

tfendif 



typedef struct{ I* Pointer AND Len */ 

ONEBYTE *ptr; 
int len; 

}PANDL; 

I* Struct and global data definitions 
*/ 



struct T_FLOAT_EPHEMERIS { 
UNSIGNED_ON EBYTE iode; 



/* Field name abbreviations were extracted from ICD-GPS-200. 



*/ 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 



toe; 

mO; 

delta_n; 

e; 

sqrta; 

omega_0; 

iO; 

w; 

omega_dot; 

i_dot; 

cue; 

cus; 

cre; 

ers; 

cic; 

cis; 

tgd; 



/* Time of ephemeris. * I 

r Mean anomaly at refernce time. */ 

r Mean motion difference from computed value. 
I* Eccentricity. *1 

I* Square root of semi-major axis. *1 

I* Longitude of ascending node of orbit plane 
at weekly epoch. */ 

/* Inclination angle at reference time. */ 

/* Argument of perigee. */ 

I* Rate of right ascension. */ 

I* Rate of inclination angle. */ 

r Amplitude of the cosine harmonic correction 
term to the argument of latitude. */ 

I* Amplitude of the sine harmonic correction 
term to the argument of latitude. */ 

r Amplitude of the cosine harmonic correction 
term to the orbit radius. *7 

r Amplitude of the sine harmonic correction 
term to the orbit radius. *1 

I* Amplitude of the cosine harmonic correction 
term to the angle of inclination. */ 

r Amplitude of the sine harmonic correction 
term to the angle of inclination. * I 

!* Estimated group delay differential. *1 



*/ 
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double toe; 
double at 2; 
double afl ; 
double afO; 
double n; 
double e_sqrt; 

double wkl ; 
double wkO; 
double sin_w; 
double cos_w; 
double a; 
double sqrta_e 

} float_ephemeris 



I* Clock data reference time. *1 

I* Polynomial coefficient (SV clock correction) 
I* Polynomial coefficient (SV clock correction) 
I* Polynomial coefficient (SV clock correction) 
I* Corrected mean motion. *1 

I* Square root of eccentricity^ subtracted 
from 1 . *1 

I* Omega_dot - WE. */ 

I* Omega_0 - WE * toe. *1 

r Sine w. */ 

/* Cosine w. * I 

I* Semi-major axis. *1 

fconst; r Square root of semi-major axis * 
eccentricity * FCONST. */ 



*1 
*1 
* / 



struct T_RECEIVER_CHANNELS { 
UN SIGNED_ON EBYTE svid; 
UNSIGNED_ON EBYTE mode; 

UN SIGN ED_ON EBYTE strength; 
UNSIGNED_ON EBYTE flags; 

}; 



struct T_GEODETIC { 

TWOBYTE degrees ; 

UNSIGNED TWOBYTE minutes ; 
double seconds ; 

}; 



struct T POS CHAN STATUS { 
UNSIGNED_ON EBYTE month; 
UNSIGNEDONEBYTE day; 
UNSIGNED_TWOBYTE year; 
UNSIGNEDONEBYTE hours; 
UNSIGNED_ON EBYTE minutes; 



double 

struct T_GEODETIC 

struct T GEODETIC 

double 

double 

double 

double 

double 

UNSlGNED_ON EBYTE 
UNSIGNEDONEBYTE 
UNSIGNED ONEBYTE 



}; 



seconds; 
latitude; 
longitude; 

datum^height; /* meters */ 
msl_height; /* meters */ 
velocity; /* m/sec */ 
heading; /* degrees *1 
current_dop; 
dop_type; 
visible_sats; 
sats_tracked; 

struct T_RECEIVER_CHANNELS channel[NUM_CHANNELS] 
UNSIGNED ONEBYTE revr status; 



#endif /* _GPSTRUCT_H */ 
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J. HEAD-C.H 



/* moxa-c.h ver 


1.00 10/11/1991 


* Definitions for MOXA Serial I/O Controller ioctrl 
*/ 


1* BAUD rate setting */ 


//define B50 


0x00 


//define B75 


0x01 


//define B1 10 


0x02 


//define B134 


0x03 


//define B1 50 


0x04 


//define B300 


0x05 


//define B600 


0x06 


//define B1200 


0x07 


//define B1800 


0x08 


//define B2400 


0x09 


//define B4800 


OxOA 


//define B7200 


OxOB 


//define B9600 


OxOC 


//define B19200 


OxOD 


//define B38400 


OxOE 


//define B57600 


OxOF 


1* MODE setting 


*1 


//define BIT_5 


0x00 1* Word length define 


//define BIT 6 


0x01 


//define BIT 7 


0x02 


//define BIT 8 


0x03 


//define STOP 1 


0x00 r Stop bits define */ 


//define STOP_2 


0x04 


//define PEVEN 


0x18 1* Parity define */ 


//define P ODD 


0x08 


//define P SPC 


0x38 


//define P MRK 


0x28 


//define P NONE 


0x00 


I 9 MODEM CONTROL setting */ 


^define C_DTR 


0x01 


^define C_RTS 


0x02 


1* MODEM LINE 


STATUS */ 


#define S_CTS 


0x01 


^define S DSR 


0x02 


#define S_RI 


0x04 


^define S CD 


0x08 



I* function declaration */ 
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long siojqueue(int); 
long sio_ifree(int); 
long siooqueue(int) 
long sio_ofree(int); 



APPENDIX F: NPS GEODETIC MARKS 



SPANAGEL GEODETIC MARKS 

WGS-84 



Maries 



SPANAGEL HALL 



2 3 1 



Level 8 



> 50.85 Deg 



Spanagel Hall runs SW to NE. The marks are on the 8 th 
level near the North-East cornor of the building. 

Marks 1 and 2 are standard survey disks, marked SPA-1 and 
SPA-2. SPA-3 is a square of thin aluminum. All marks 
are attached to the roof under a rail. 

Poles are mounted on the r6il above the marks. The 
height to the marks are 11.83 m ( WGS-84 elipsoid ). 



Mark Pole 



SPA-1 


36 35.70375 


121 52.47777 


11.83 


14.08 


SPA-2 


36 35.70180 


121 52.48074 


11.83 


14.10 


SPA-3 


36 35.70277 


121 52.47926 


11.83 


14.09 
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Table IX WGS 84 SPANAGAL HALL EIGHTH FLOOR GEODETIC MARKS. 



Loca- 

tion 


Item 


degrees 


degrees 

minutes 


degrees 

minutes 

seconds 


meters 

(above 

geoid) 


Spa- 1 


Lat 


36.59506253 

°N 


36° 35.70375’ 


36° 35’ 
42.2251” 




2321 


Lon 


121.87462956 

°W 


121° 

52.47777’ 


121° 52’ 
28.6664” 






Alt 








11.83 


Spa-2 


Lat 


36.59503006 

°N 


36° 35.70180’ 


36° 35’ 
42.1082” 




2322 


Lon 


121.87467903 

°W 


121° 

52.48074’ 


121° 52’ 
28.8445” 






Alt 








11.83 


Spa-3 


Lat 


36.59504622 

°N 


36° 36.70277’ 


36° 35’ 
42.1664” 




2323 


Lon 


121.87465439 

°W 


121° 

52.47926’ 


121° 52’ 
28.7558” 






Alt 








11.83 



140 



LIST OF REFERENCES 



Advantech, PCLS-802 PC-ComLIB Serial Communication Programming 
Library, Revision 3.0, December 1992. 

Anderson, John D., Jr., Introduction to Flight, New York, 
McGraw-Hill Book Company, 1978. 

Applied Research Laboratories, GPS Seminar, presented to the 
Naval Security Group, University of Texas, Austin, 26-28 
January 1993. 

Beck, G. E. , Navigation Systems: A Survey of Modern Electronic 
Aids, Van Nostrand Reinhold, 1971. 

Bowditch, Nathaniel, American Practical Navigator: An Epitome 
of Navigation, Defense Mapping Agency Hydrographic/Topographic 
Center, Vol I, 1984. 

Carey, Susan, Sextants in Space Can Change the World, Wall 
Street Journal, v. CXXIX, no 13, p. Bl. 

Chan, Ling, "GPS World Receiver Survey," GPS World, v. 4, no. 
1, pp. 52-65, January 1993. 

Dunlap, G. D., and Shufeldt, H. H. , Dutton's Navigation and 
Piloting, 12th ed., United States Naval Institute, 1969. 

Hobbs, Richard R. , Marine Navigation 2: Celestial and 

Electronic, Naval Institute Press, 1983. 

Kayton, Myron, Navigation : Land, Sea, Air and Space, New York, 
IEEE Press, 1990. 

Kremer, G. T. , and others, The Effect of Selective 
Availability on Differential GPS Corrections , Navigation, 
Journal of the ION, v. 37, No. 1, Spring 1990. 

Leick, Alfred, GPS Satellite Surveying , New York, John Wiley 
& Sons, 1990. 

Lin, Ching-Fang, Modern Navigation, Guidance, and Control 
Processing, pp. 215-231, 445-464, Englewood Cliffs, N.J., 

Prentice-Hall, Inc., 1991. 

Logsdon, Tom, The Navstar Global Positioning System, Van 
Nostrand Reinhold, 1992. 



141 



Motorola, GPS Receiver Technical Reference Manual, Revision 
4.0, 1992. 

U.S. Naval Oceanographic Office, H.O. Pubs 117A and 117B Radio 
Navigation Aids, Government Printing Office, Washington, DC, 
1972 . 

Rockwell International Corporation, Interface Control Document 
ICD-GPS—200 (Revision B) , 1987. 

Tattersf ield, Donald, Orbits for Amateurs with a 
Microcomputer, New York, John Wiley & Sons, 1984. 

Wells, D.E., and others, Guide to GPS Positioning, 2nd ed. , 
Fredericton, N.B., Canada, Canadian GPS Associates, 1987. 

Wilkes, Owen and Gleditsch, Nils Petter, Loran-C and Omega: A 
Study of the Military Importance of Radio Navigation Aids, 
Oslo, Norwegian University Press, 1987. 



142 



INITIAL DISTRIBUTION LIST 



No. Copies 

1. Defense Technical Information Center 2 

Cameron Station 

Alexandria VA 22304-6145 

2. Library, Code 52 2 

Naval Postgraduate School 

Monterey CA 93943-5101 

3. Professor Michael Shields, Code EC/SL 4 

Department of Electrical and Computer Engineering 
Naval Postgraduate School 

Monterey, California 93943-5121 

4. Professor Isaac Kaminer, Code AA/KA 2 

Department of Aeronautics and Astronautics 

Naval Postgraduate School 
Monterey, California 93943-5121 

5. Professor Rick Howard, Code AA/HB 1 

Department of Aeronautics and Astronautics 

Naval Postgraduate School 
Monterey, California 93943-5106 

6. Professor Robert McGhee, Code CS/MZ 1 

Department of Computer Science 

Naval Postgraduate School 
Monterey, California 93943-55118 

7 . Lieutenant Commander Eric Twite 2 

32 Bahama Bend 

Coronado, California 92118 



143 



DUDLEY juMOa lsBRARY 
'•f'AV' ‘ 

IVfGNl W. .. 



'^n’UfVV 

.*'1 rfww'i 

i-1 



